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ABSTRACT 


Cebeci’s viscous inviscid interaction program was applied to the analysis of 
steady two dimensional incompressible flow past a NACA 65-213 airfoil at zero angle 
of attack at a Revnolds number of 240,000. Predicted boundary laver characteristics 
Were found to be quite sensitive to the choice of boundarv laver transition begin and 
length. Good agreement with the experimental results of Hoheisel er al could be 


obtained by proper choice of both transition begin and length. 
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I. INTRODUCTION 


The prediction of the stall characteristics is of major importance in aeronautical 
engineering to determine the operating limuts of an aircraft. Therefore, the 
development of reliable and accurate numerical methods for predicting separated flow 
regions is one of the most challenging problems in Computational Fluid Dynamics 
(CFD). 

In this thesis we lumit ourselves to the problem of incompressible two- 
dimensional airfoil flows. Two basic methods are available to compute viscous flows 
Which include regions of flow separation. The first approach is based upon a solution 
of the full Navier-Stokes equations (or some approximate form, such as the parabolizéd 
Navier-Stokes equations). This approach has the disadvantage of being very expensive 
and time consuming. The second approach 1s based upon the so-called Viscous-Inviscid 
[nteraction Method. The outer flow is computed using the inviscid flow equations. The 
inner flow (close to the airfoil) is obtained: from a numerical solution of the Prandtls 
boundary laver equation. However, in contrast to the well-known classical boundary 
layer computations the pressure cannot be prescribed a priori, but must be found 
iterativelv (4e., by viscous-inviscid interaction). This approach has the advantage of 
being much faster and more efficient than the Navier-Stokes solutions. | 

The approach chosen in this ‘thesis is based upon the viscous-inviscid interaction 
method developed by T. Cebeci and collaborators at the Douglas Aircraft Company. In 
chapter 2 the fundamental equations are sunimarized . Chapter 3 1s devoted to a 
discussion of the inviscid flow method, especially the so-called Panel Alethod first 
introduced [Ref. 4] by Hess and Smith. In chapter 4 the direct and interactive 
boundary laver methods are discussed, followed by a brief explanation of the 
turbulence model used. Finally, the computer program is explained and computed 
results are presented for the NACA 65-213 airfoil and compared with detailed 


measurements by Hoheisel et al. [Ref. 14]. 
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~ TI. FUNDAMENTAL EQUATIONS 


This chapter presents the equations used in our analysis, which involve the 

development of: 

[. Conservation of mass (Continuity equation) 

2. Conservation of momentum (Newton's Second Law) 

>. Inviscid rlow equations 

4. Boundarv laver equation 

Pee routcice niodel 
With these cyuations we can predict the behavior of a body moving through the fluid. 
Mme@erecase We ire dealing with an airfor in two dimensional. steady. inviscid and 


MISeOUS (ON. 


A. CONSERVATION OF MASS 

Let us apply the principle of conservation of mass to a small volume of space 
Geaee re cate {uid can move ireely. For convenience, we shall use a cartesian 
Seememiatte ssien ix. .c). Furtherniore. in the interest of simplicity, we shall treat a 
->-D flow, that is. one in which there is no flow along the = -anis. Fiow patterns are the 
same for anv .x-) sane. avs Indicated in thessketch of Pieure 2.1. the component of the 
fluid velocity in the x direction will be designated by uw, and in the v direction bv v. The. 
nec outflow of mass through the surface surrounding the volume must be equal to the 
Meeicedseso, Imass wittin the Volume. The mass tiow rate through the surface bounding 
Mierelemeniisecqual toe the product of the density, the velocity component normal to 
ieeserrtee, and area of that surface. 
pmimst-oncer Pavior sevies expansion is used to evaluate the flow properties at the faces 
of the element [Ref. 3] since the properties are a function of position. Consider the 


flow out of the volume as positive, then the net outflow of mass per-unit ume 1s the 
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Figure 2.1 Sketch wdlustrating the velocity 
and the densitv for mass flow balance 
through a fixed volume in 2-D 


which must equal the rate at which the mass contained within the element decreases 


Op 9@ a 
at gel Pl > sor) = 0 (2.3) 


femnec-of form. eguacion 2.5 is 
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PieetccmemenoreSsure Variations that occur in relatively low speed flow are sufficiently 
pee me censi.» (Ss essentially constant. For these incomopressibe flows, the continuity 


equation becomes 


a 


Or Oy 


or. Ip vector form V-} =0 


Ee eee nse ny ATION OF NTONMENTUMI 

[ee craton of tie Conservation of finear momentum is obtained by appiving 
Re eee ree, >| where the net force acting on a fluid particle ts equal to the 
time race of change of the linear momentum of the fluid particle. As the fluid moves in 
Seeeee sishace dnd volume may Change, but its mass is conserved. Thus, using a 
CoOnemmes Sistem that is neither accelerating nor rotating. called an inertial coordinate 


So stequerre mat write 


DY 


= m— (2e7) 
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The velocity } of fluid particle is. in general. an explicit function of tume t as well as 
of its position xy. Furthermore, the position coordinates x, of the fluid particle are 
themselves a function of time. Since the time differentiation of equation 2.7 follows a 
given particle in its mouon, the derivative is trequentin= ternicu em pan elcmmon 


Substantial derivative On mc iieemmew ew nieagit ny.) mr (a 
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Thus, the substantial derivauve is the sum of the local, time desevuenmenemeas 
that occur at a point in the flow field and of the convection in space. “whem the loaae 
time dependent changes are zero, GJ’ c:=0 , such flows are known as steady-state 
flows. The principal forces that act en the body are these whien act directly on ine 
mass of the fluid element, the body forces , and those which act on its surface. the 
pressure forces and shear forces known as surface forces. The stress svstem acting on an 
element of the surface 1s illustrated in Figure 2-2. The stress components ft acting on 
the small cube are assigned subscripts. The first subscript indicates the direction of the 
normal direction to the Surface ones Miceli the stress caer sean mG Gmce ae Meme Miosemian 
Insicates the direction in which the stress acts. Thus. t,,. denotes a stress acting in the 
v direction on the surface whose normal points in the \ direction, The properties of 
most fluids have no preferred direction in space; that 1s. fluids are isotropic. 

Again, the stresses on the fluid element can be obtained from a Tavlor series 


expansion. In general, the variolis stressesvenineeeinemiepo@ini slo POImtNy Tiniss ties, 
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Figure 2.2 Stresses acting on a 2-D element of fluid. 


produce net forces on the flumd particle. whieh cause it to accelerate.) ie 1 oreeseenme 
on each surface are obtarned By Wiking into deeem@nt the variaiemsolesicsom ain 
position, by using the center of the element as a reference point. 

To simplify the illustration of the force balance on the fluid particle we shall 
again consider a 2-D flow, as imdicated in Figure 2-2.eeiine nesultantsiGr ce sige 


direction. for one uni lemenin mas 
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Where f is the body force per-unit mass in the x direction. The most common body 
force for the flow fields is that of eravit.: Equauen 2310) is the leitiancesidemen 
equation 2.7. Por the pelt hand side. vomibimtesmidsseteriiy iti ec Ualioneen snimnesss 


direction 
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From equation 2.10 and) 2:11) substuture 1ntoreq inc meer divided bi sername 
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Next. we need a relation between the stresses and the motion of the fluid. Fora 
fluid at rest or for inviscid fluid motion, there 1s no shearing stress and the normal 
Stresses in the naturcof aupressupes or IMG pat (ele men emo une cc mioerciited [Omics are 


of strain bv a physical law based on the following assumptions: 


L6 


1. Stress components mav be expressed as a linear function of the components of 
the rate of the strain. The friction law for 1-D flow of a Newtonian fluid is a 
special case of this linear stress;rate of strain relation, ie., tT = p ( Gu dy ), 
where [If 1s the fluid viscosity. 

2. The relation between the stress components and strain rate components must 
be invariant to a coordinate transformation consisting of either a rotation or a 
murror reflection of axes, since a physical law cannot depend upon the choice of 
the coordinate system. 


When all velocity gradients are Zero (1.e., the shear stress, vanishes), the stress 


(ss 


components must reduce to the hydrostatic pressure p. 


For a fluid that satisfies these criteria. in two dimensional flow 
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with the appropriate expressions for the surface stresses, substitute into equation 2.12 


and 2.13, we have 
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These general differential equations for the conservation of linear momentum are 
known as the Nuvier - Stokes equations. When we are dealing with incompressible and 


2-D flow, then fron equation 2.6, V.F" = @, and the body forces are neglected, and t., 


+ t., = - 2p. Therefore equation 2.17 and 2.18 can be written as 


u a0 +y Ou on 1 Op O* meen 
se a rp} — 9 
Or Oy Ot pdr | dpe dy? (2.19) 
Ov Ov Or 1dp 02» p2p 
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These are the well known Navier - Stokes equation for 2-D, incompressible, viscous 


flow. 
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C. INVISCID FLOW EQUATION 

inviscid flow represents an ideal flow, where the eliects of viscosinn sane rzenomeia 
reality, this 1s not true because every medium has viscosity, eventhough it may be very 
small. | 

Why is the inviscid flow important in dealing with a body moving in the viscous 
fluid ? 

In 1904 L.Prandtl came up with the answer. For high Revnolds Number on a 
bodv moving in a viscous fluid, two regions can be distinguished. The effect of 
viscosity can be neglected outside a ver thin region mear the body“ hiciwiss< tee 
boundary layer. Inside the boundary laver the viscous effects are important (further 
discussion in viscous flow section). This is the reason why the inviscid flow remains 
Important in the computation of fluid dvnanucs, eventhough it represents an ideal case. 

|. Potenial Flow 

Since the flow upstream of the bodv is uniform then it is also irrotational ( in 


inviscid flow). 


Consider the following vector identity , if @ is a scalar function, then 
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1.e., the curl of the Gradient of a scalar function is identically zero. Comparing equation 
See edie oe. Wie See that 


— 
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Equauon 2.25 states that for an irrotational flow there exists a scalar function @ such 
that the velocity is given by the gradient of @. We denote @ as the velocity potential. © 
is a function of the spatial coordinates. 1.e..@ = (xv). And from the definition of the 


gradient in cartesian coordinates . equation 2.23 can be written as 


8 
Le and c= Sy eteua 
Os (2.24) 


Tnus, @ has the property that its parual derivauye in any direction is the 
velocity component in that direction. It follows that the existence of @ is the sole 
criterion for irrotationality. The usefulness of the elOctverotential imeilous tof 
Podemedisienilicance derives from: the circumstance that. for a bodv in relative motion 
in an originally irrotational flow. the circulation vanishes around any contour that does 
not include the bodv or does not intersect the boundarv laver or the wake. therefore, a 
velocity potential can be found to describe the flow evervwhere outside the boundarv 
ieeeror tie wake. When a [low field is irrotational, hence allowing the velocity 
potential to be defined. there 1s a tremendous simplification. Instead of dealing with the 
velocity components (say, uv and w) as the unknowns, hence reyuiring three equations 
[Omeemeestignenowns., we Can imsteau deal with the velocity potential as one unknown, 
therefore, requiring the solution of only one equation for the flow field and the velocity 


components can be obtained from equation 2.24. Because irrotational flows can be 


described by the velocity potential @ such flows are called potential flows. 
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-. Governing Equation for irrotational. incompressible flow : Laplace Equation. 
We have seen from equation 2.6, that the consermmation of niassm@for an 
incompressible flowtakes the form V.i’ = 0. In addition. for irrotational flow we have 
Seen Ih eqiuale imme m7oam y= Vo. Therefore. for a flow that is both incompressible and 


also irrotauional. equations 2-OvanGdls-. ca bercomeimcustem: (olus 


rr u — ee 
Ss = or apt ye 0 ) 


Equation 2.25 is called Laplace’s Equation . one of the most extensively studied 


equations in mathematical physics. ° 

ete tlie Laplace's equation 18 a second onder linear partial Cilterenum 
equation. The fact that it 1s dear is particularly important. because the superpositien 
of anv particular solution of a linear differential equation is also a solution of the 
equation. For example, if @,. @,. @....., represent n separate solutions of equation 
2.25. then the sum @ = @,  @, + .... + @, is also a solution of equation 2.25. 

Since irrotational. incompressible flow 1s governed by Laplace's equation and 
Laplace's equation 1s dinear, we conclude that a complicated) Mow sattenimmonecde 
irrotauonal, incompressible flow can be synthesized by superposition of elementary 
flows which are also irrotauonal and incompressible. The singularity (or panel) 


methods presented im the next chapter are pbaseduonmnume ca. 


D. BOUNDARY LAYER EQUATION 

Up to this point. we have been dealing with tien outside the boundanm slavem 
where viscous effects remain small. In the region within the boundary laver, velocity 
gradients are high even with vers smiall viscosity munenciore. 16 Dec Olicss.eis inp encaan 
to deal with a real fluid. Before the boundar viayercan be analy zed iittaers wierd. 2 
to know what governing equations can be used in the practical analvsis. The objective 
IS tO predict viscous flows by means of the boundarv laver method, instead of solving 
the complete Navier-Stokes equations. 

From the previous derivations, equdieis 2.5 meme) (numer Omdnemiscu licence sun 
Orderto simplify these equations we have to make some assumptions: two- 


dimensional, steady, constant fluid properties, and no body forces. Another important 


assumption 1s that the boundarv laver thickness is very small compared to the length of 
the body (airfoil in this case). With these assumptions, L. Prandtl in 1904 introduced 
the “order of magnitude” estimate into equation 2.5, 2.19 and 2.20. In order to do this, 
all linear dimensions will be referenced to the characteristic length / and all velocities 
Pee eccierenced foc. theretone (and € can be said to have order of magnitude of 
one ( Written as O(1) ) and with the above assumptions, vy will correspond to the 


boundary laver thickness (6). Then the continuity equation 2.5 can be written as 
U U6 
(7) (3) a () = oF) 


Because 6 is Verv small. from the above relation we can imply that v, the velocity in 
the v direction must be very small. Therefore. the contuunuitv equation sull holds in the 


boundary laver. From the Navier-Stokes cquation in the x direction (eqn 2.19) 
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[ntroducing the order of magnitude in the above equation. we have 
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If all above terms multiplied by ( then 


O) (1) (za) 


( 1 [2/52 
Ul/v | Ui/v 
[f the Revnolds Number is high, the fourth term 
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For the v direction ( equation 2.20 ) in terms of order of magnitude we obtain 
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The only term left is the’pressure term, because dimommervemms are oO! fischer Once 


ay : (2.27) 


This 1s very important because it tells us thatetne mpessimmenacross the Doundan wave: 
remains constant. The triplet of equations (925652 20 and 2.2) )sanueticupevigucms 
condition of zero normal and tangential velocity are known as the Prandtl’s boundary 


layer equations. 


FEF. TURBULENT FLOW EQUATION | 

Since the continuity equation (2.5) and the Navier-Stokes equation (2.19) make 
no assumptions regarding the tvpe of flow, thev are instantaneously valid in both the 
laminar and twrbulent flow regimes. However, it ts too difficult to deal with 


instantaneous properties in turbulent flow. Therefore we introduce the tme-averaged 


tJ 
tJ 


properties. In our notation, prime denotes the fluctuation quantity, and bar denotes 
the mean value. 
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Introducing these properties into the continuity equation and averaging, we have 
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Garrvine out the integration term by term (details are in {Ref. 17] }. vields 
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Therefore the continuity equation fér turbulent flow becomes: 
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oe an (2.28) 
A simular procedure can be applied to the Navier-Stokes Equation 
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The last two terms correspond to the normal and shear stress terms respectively, which 
we call the Reynolds Stresses or Turbulent Stresses. 
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If. PANEL METHOD 


A. INTRODUCTION 

The thin airfoil theorv is just what it says, it applies only to a thin airfoil at small 
angle of attack. It is not much used these davs for the analvsis or design of single 
element airfoils. It does give fairly good results for airfoils of 12 °o thickness or less. 
On the other hand, the determination of the aerodynamic characteristics of thick. 
highly cambered, slotted surfaces, with single or multiple flaps and mutual interference 
effects among wings, fuselage. nacelles, and so forth, requires, in general, the use of 
numierical methods. ; | 

In the 1960’s Hess and Smith at McDonnell Douglas introduced a method. the 
so called Punel \ethod (Ref. 4], as a numerical approach for 2-D flows which can be 
extended to 3-D potential flow problems. Such methods are called panel methods 
because the body surface is approximated bv a collection of panels . There are a 
number of wavs to set up the panel method. io begin with, there are choices evened. 
(0 the tipesor singularity used, sources, doublets, vortices or a combination of source 
and vortex distributions. 

Panel methods as a numerical approach for predicting forces and moments acting 
on the body gave good agreement with the reliable published data. The applicauon of 
the panel method requires that the problem can be formulated such that 

1. the body can be represented bv a closed polygon of a finite number of elements, 


called panels connected by nodes. 


tw 


the flow tangency condition is satisfied in the middle of the panels (control 
points) to avoid the inaccuracies of thin airfoil theorv. 

3. the singularity distribution of each element is approximated by some kind of 
analytical function. Also, the singularities should be distributed on the body 
surface rather than on the chord line or anv other line within the body. 

Before we discuss further details, it is necessary to know about the basic formulation 
for source and Vortex distribution as a singularitv parameter. 
|. Single Source and Source Panel 

A flowfield where all streamlines are straight lines emanating from a source 


point, is called a source flow. On the other hand, when the flow direction is inward, it 


1s called sink flow. V.I’ = @ at every point except at the origin where V.I° becomes 


infinite. The source flow is irrotational at every point. Its velocity potential is defined 
as a 


single source Wea: a 


where .\ is defined as the source strength and r is the distance from the considered 
point (xg) to the source. When we are dealing with non-lifting flows over a bodv, we 
can superimpose elementary source fiows in order to obtain a complete solution. This 


method 1s called source pdnel method. 
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Meme! is the panel lencth. 
-. Single Vortex and Vortex Panel 
The flow where all streamlines are concentric circles about a given point is 
velined as single vortex flow. The velocity along anv given circular streamline is 
C@retadv cul varies (rom Streamline to streamline. Its velocity potential can be written 
as 


r 
single vortex cA fe me ee a G Bo) 
ym ye . 
where @= arctan] ———/ | with (++.y,-) as the center of vortex 
mS 


etmis wmae@ine a straight dine perpendicular to the page (Figure 3.1). This line 1s a 
straight vortex filament of strength T’ and the flow induced in anv planes perpendicular 
to the vortex of strength [. 1.e., the flows in the plane L to the vortex filament at o 
Sere “ane identical ‘0 e€ach other and are wentica!l to the flow induced by a point 


vortex of strength [. 
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Figure 3.1 Vortex Sheet Representation. 
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Refer to Figure 3.1, imagine an infinite number of straight vortex filaments 
side by side, where the strength of each filament is infinitesimally small. Define y = 
¥(s) as the strength of the vortex sheet per unit length along s, then the strength of the 
infinitesimal portion ds of the sheet is yds and the small section of the vortex sheet of 


strength yds induces an infinitesimally small velocity dv at point P(x,z 


ds 
di = "“75r otto er ‘ 
“0 (3.4) 


It is sometimes more convenient to use the velocity potential @, and the increment in 


velocity potential dw induced at P(x,z) bv elemental vortex yds 


For the entire vortex sheer, froma to b 
y 
cian — =| O~ds | 2 
tiers (3.6) 
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Equation 3.5 is useful in classical thin airfoil theory and equation 3.6 is important for 


we 


the numerical vortex panel method. 


B. SOURCE AND VORTEX DISTRIBUTION 
Having introduced the basic idea of the panel method, we can use these 
singularities either alone or in combination. This method ts due to Hess and Smith. 


Thus, the potential mav be decomposed in a manner such that 


t 
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with @, being the potential of the uniform onset flow, and . and @, the potentials due 


to source and vortex distributions respectively, hence 
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| per ae (3.8) 


ae 
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in which the integrations are to be performed over the body surface. Because of the 
superposition principle, this @ automatically satisfies Laplace's Equation ( see equation 
2.25) and the boundary condition at infinity. It will be the solution we seek. if o(s) and 
“(s) are deternuned so as to meet the boundary condition of flow tangency and the 
Autta condition (to be discussed in the next section). 

Hess and Smith assumed the vortex strength to be constant over the body 
surface and the source strength must Vary Over ihe Surlace Sincesthe KUm@scondiien 
involves oniv the trailing edeeé, the vortex Strength cam Ge represented ay a sinaie 
number. Thus, if one distributes on or within the body surface. vortices whose net 
strength is the correct circulation, the problem is solved if sources can be distributed 
over the bodv surfacevsa as to miake the toral velocim Wield (compnmsed of the oases 
flow and the velocity fields due to sources and vortices} tangent to the body suriace 
regardless of how the vortices are distributed: “Homeversthe intceralsmecatme: seams 
and 3.9 are hard to evaluate. even for simple forms of the source and vortex strength. 
unless the surface on which the sources and vortices are distributed, Is a straight line. 
Thus, we select a certain number of points on the body contour, called nodes, and 
connect the nodes with straight lines, which become the panels of the method (see 
Figure 3.2 ) | | 
We then distribute the sources and vortices on the straight line panels. so that the 


potential Civen DWeaquation si) Cansace. mi comeas 


g; Bf 
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[n most cases, equation 3.10 still allows an exact solution of the flow problem. The 
exceptional cases are those in which the sources and vortices must be distributed 
exactly on the boudv surface ; to be™mianicmancaly  paeccise mgm vic meu cmp otemna| 
cannot be continued anlytically across the body surface. By increasing the panel 
density, the body shape can be better approximated. This is the only major 
approximatuon of the panel method, one that becomes more accurate as the number of 


panels increases. 


Figure 3.2 Discretization: a), Actual airfoil 
). Airfoil after discretization. 


29 


— 
Pat | 
ah 





(2, y;] 


(xo eee 


ds ; 
0; 
—_ —— / —— —— —s —s 
/ Be 
Ean 
/ 
et 


is 


—} 


Figure 3.3 Representation of i-th and j-th panels 
in the Panel Vlethod. 


To implement this method, we need a nomenclature. Let the i-rh panel be 


defined as the one between the i-rh and (+ /)th nodes, and its inclination to the x axis 
be @.. Therefore the relation between midpoints and boundary points 
fey a5 yey 


Z 
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and the computation of the angle @ 
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Pe@mretne 2eOmetry in Figure 3.5 a relation for angle and distance between two panels 


can be obtained 
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1. Flow Tangency Condition 
Tne flow tangency condition in the case of no blowing or suction, 1s satisfied 


at the middle of each panel sometimes called rhe no penetration condition. To obtain 
MOM resoeem to the mormal vector for each 
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_ 


ealsecemeaition, differentiate equation 
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ee The influence coefficient or normal velocity (normalized with Vo) at the 


midpoint of the i-ch panel due to a source distribution at the j-ch panel. 


BB;,; --- The influence coefficient or normal velocity (normalized with Vj) at the 
midpoint of the i-ch panel due to a vortex distribution at the j-ch panel. 


z,,yj-- Coordinates of panel midpoint of i-th panel 
Coordinates of panel midpoint of j-ch panel 


Tey ite Se, an length of j-ch panel 
2. Concept of the Influence Coefficient 
Equation 3.12 can be evaluated by integration. To obtain this, we have to 
relate the coordinate of a point in the j-th panel with the known coordinates of the 


boundary points and panel angles. From the geometry in Figure 3.2, one obtains 


a cos 3; (3.13) 

On; 

Orr; 
but, cos B. = - sin 8. and sin B. = cos @, , therefore equation 3.13 and 3.14 become 
gan = — sin 0; (3.15) 

Onn; 
oy = cos); (3.16) 
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On the other hand, 
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And the influence coefficients can be obtained in terms of 
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Adi = 5-[5CF - DG] (3.19) 
_ Ml 
BBi; = 5-|5DF- CG (3.20) 


where, 
AS (ey — -\,) c089, —(y, — 1 ,) sin ?, 
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5. Computation of Total Disturbance Velocity V 


The computation of the total disturbance velocity in x and y direction can be obtained 
Gree crenilatingsin the x and ¥ direction. 
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where 
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caedey The horizontal velocity component at midpoint of i-rh panel due to a unit 
source distribution at the j-rh panel. 
Ad’. The vertical velocity component at midpoint of i-rh panel due to a unit 


source distribution at the j-ch panel. 


BB, .. The horizontal velocity component at midpoint of i-rh panel due to a unit 


vortex distribution at the j-rh panel. 


Bee The vertical velocity component at midpoint of i-ch panel due to a unit 
vortex distribution at the j-rh panel. 
These influence coefficients can be obtained bv integrating equation 3.12 and 


introducing equations 3.15 thru 3.18 into equation 3.12 


BE =| ; cos8,F + sin d;G} (3.24) 
BBY = ~ = sin #,F ~ cos 8, G} | (3.25) 
ay = ee (3.26) 
BB! = ~ Ady, (37253) 
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4. Kutta Condition 

As in all problems concerning airfoils in inviscid flows, an auxiliary condition 
needs to be invoked to ensure that a unique solution is obtained. This condition, 
known as Autta Condition. relate to assumptions about the flow characteristics at. or at 
least in the neighbourhood of, the trailing edge. When the surface velocities are made 
equal at the nudpoints of the trailing edge elements then, by the Bernoulli Equation, the 
pressures at these points are also equal. so the Kutta condition can be represented as 
the condition of zero loading in the region of the trarling edge’ (Ref. 1] which is 

physically realistic. Therefore the Kutta condition can be summerized as follows: 
ipeetor a civenwvangle of attack. the value of circulation F around the airfoil is such 


that the flow leaves the trailing edge smoothlv. 


ty 


If the trailing edge angle ts finite, then the trailing edge is a stagnation point. 


If the trailing edge angle ts cusped. then the velocities leaving the top and the 


os 


bottom surfaces at the trailing edge are finite and, equal in magnitude and 
direction. 
These are the basic principles of the Kutta condition. In the actual comiputation we are 
using in the code. there is no restriction whether the trailing edge angle is cusped or 
not as mentioned in item 3, but.rather, two tangential velocities in the last control 
points assume have the same magnitude but in opposite direction. 

In the numerical solutions. the actual tratling edge is not a stagnation point. 
Furthermore, it is found that the velocities at the nudpoint of the trailing edge elements 
differ significantly from stagnation values; they are more likely to be closer to the free 
Stream velocities. 

5. Determination of the vortex strength 
The Kutta condition is used to determine the vortex strength ¥. From 


equation 3-12 , set ¥ = 1.0 and obtain the expression of oF in terms of the influence 


coefficient 
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Nv Vv 
Le Adi oy = | Ee BB,,|-|sin(4 - a a 


j=l 721 


This equation canbe solved by using Gaussian elimination. The result can be 
substituted into equation 3.19 and 3.20. By letting ¢=/ and :=.N for the trailing edge 
panel, we can compute horizontal and vertical velocities for i-th and .V-th panel as ae 
Mop Vey and V.\. Introducing these values into the Kutta condition, produces a 
quadratic equation in ¥y. This method is not the best wav to satisfv the Kutta condition 
in steady flow calculation, but by using this method the code can be extended to the 


unsteady case. 
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a, and B. are the result obtained from equation 3.28. From equation 3.29, there are 
two values of vortex strength 7, but only one value is used in the further calculation. 
These values relaté to the angle of attack. If the angle of attack is positive, use the 
negative value of y. If angle of attack is negative use the positive value. 
6. Determination of the source strength 
Once vortex strength has been determined using equation 3.12. the source 


Strength 6 can be obtained from 
[o|= s{aj] + [3] (3.30) 
7. Calculation of ‘on bodv’ velocities 


The velocity at midpoint of the i-ch panel can be obtained by a spatial 


derivative of the velocity potential in the tangential direction. 


° 0 27 94 
\; = a7 el ti vill (3.31) 
SS « | 
i= LAT oj; +7 L BT; + oslo; - a) (522) 
: j=l j=l 
ee | Pe 
where iAT,; = --(-CG— =DF| (rss) 
1 f ‘ 1 a4 
Pap =| Dee C | (3.34) 
Mem i= j Ale = 00 (3.35) 
It 
Coc iaas (3.36) 


AT,, .1s the influence coefficient of tangential direction at the midpoint of the i-th 
panel due to a unit source at j-th panel, and B7,, is the influence coefficient of 


tangential direction at the midpoint of the i-th panel due to a unit vortex at j-th panel. 
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a. Calculation of Pressure Coefficients 
Once the velocity on the body has been calculated, we can easily obtain the 


pressure coefficientS-by the Bernoulli equation 


= = wwowle (3.37) 
sA\ iy 


Where V. is the dimensionless velocity (normalized by Vj) at the midpoint of the :-ch 
panel. Furthermore, as soon as the pressure coefficients are known other coefficients 
such as Lift Coeflicient(CL), Drag Coefficient(CD) and Moment Coefficient about the 


leading edge (Cf) can be obtained bv using pressure integration along the bodv 
contour which is approximated bv a closed polygon. 
= (eos) 
C; a zs Cree a tall ee! 
NV 
Cy = os — 7,,{( Xai - Ail (So) 
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C. LINEARLY VARYING VORTEX DISTRIBUTION 
The following method is only one variation of the use of the panel method 


[Ref. 7] ; it involves representation of the airfoil by a closed polvgon of the vortex 


panels. [he vortex-panel method introduced here has the feature that the circulation 
density on each panel varies linearly from one corner to the other and is continuous 
across the corner as indicated in Figure 3.4. The airfoil and wing problems can be 
solved by means of a vortex-panel distribution alone, but calculation of fuselage and 
nacelle characteristics and their interference flows dictate the use of source and. 
possioly, doudiet as well as vortex-panel!s. For steady flows, it has been found that th 
use of piecewise constant. discontinuous distributions of vorticity can lead to 
inaccuracies such as oscillating values of the vorticity on successive panels. The use of 
a nearly varving. continuous distribution of vorticity eliminates this proolem. 

In the presence of unuorm flow Va at an angle of attack @ and m vortex panels. 


the veiocity potential at the i-ch-control point (x...) 1s defined as 





= Sus teen ( 3.42} 
Sf als) Ln = 9) 

fle.) =lele,cosa-ysinal- | — seen | Ai ds, {3.44) 
| j=l ee ee, = el 

a . _ (3.45) 


where w= 7-4 y si 
J 


is che vortex distmbution which is linear along the panel and continuous across the 
Jommdary points. . 

Tne panels, \ in number, are assumed planar and are named in the clockwise direction, 
Starting from the trailing edge, and boundarv points selected on the surface of the 
a Oimaresthe intersections of continuous vortexpanels. The (.\+/) values of 4; at 
boundary points are the unknowns to be determined numerically. The condition that 
the airfoil be a streamline is met approximately at control points. The boundarv 
condition requires that the velocity in the direction outward normal vector n. be 


vanishing at the i-th control point, such that 


ease 0 a (0) on eo ( 
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Figure 3.4 Replacement of an airfoil bv vortex panels 
of linearly varving vortex strength. 
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Carrying out the calculation of equation 3.44 and normalizing the vortex strengths by 2 
m V9 we have 


— 


ty a Pa 
t 3; 6] t~ J) 3.4 
| lg eee) 7 On; arcrau| en |, = 
2 


i) 


and from differentiation and integration we get 


Vv 


(ON, 4) + ON. 41] = sin Gi - a) (3.48) 
J=1 


Where, Y = Y'2mVpy Is a dimensionless vortex strength. The coefficients in the 


Sanemeneses are 
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4= -—(2; — -Y,)cos4; — (y; — Y,) sin 4; 
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The expressions in the parentheses on the left side of equation 3.48 represent the 


normal velocity at the i-¢h control point induced by the linear distribution of vortices 


4] 


on the j-th panel, called influence coefficient of normal velocity (except for i = m+ /). 
When / = /, the coefficients have simplified values 

ONT. = =] 

CN 


oe 
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which describe the self-induced normal velocity at the i-th control point. 


lL. Kutta Condition 
In order to apply the Autta condition in the theoretical analvsis, we need to be 


more precise about the nature of the flow at the trailing edge. It was mentioned before 
that the trailing edge can have a finite angle or it can be cusped. The statement of the 
Kutta condition in terms of the vortex sheet is y (TE) = V_- V, where Vy, and V, are 
the tangential velocities at the upper and lower side of the trailing edge. However. for 
the finite-angle trailing edge, the velocities at the upper and lower surface have same 


magnitude. Thus the Kutta condition becomes 7 (7E) = 0, 


° 


Po Pols — 
Piney =O (3.49) 


For computer programming, equation 3.48 can be arranged such that the: equation is : 


easv to formulate in matrix form ; 


iV 


0, = RS, ia) oc ee (3.50) 


J=1 


After combining equation 3.49 and 3.50 they are sufficient to solve for the (V+ /) 


unknown 1; values, in which the influence coefficients can be classified as 


108 1 ee Coe 
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From the above equation, it can be expressed in matrix form 
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To obtain the unknown 1; values, the above expression can be solved by using the 
Gaussian Elimination method or anv other linear equation algorithm. 
-. Calculation of the ’on body’ velocities 
The unknown vortex strength having been determined, we can proceed to 
compute the velocities and pressures at every control point. Recall the no penetration 
condition at the control point. Hence the disturbance velocity induced at control points 
is only the tangential velocitv. The velocity can be obtained bv differentiating equation 


3.44 in the direction of the tangential vector on the i-th panel. hence 


3 pee 
Ve = = a is i Voila 
ae An ( yi) 
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in which the coefficients in the parentheses are defined as 
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The expression in the parentheses following the summation symbol has the physical 
meaning of the tangential velocity at the i-th control point induced bv the vortices 
distributed in the j-ch panel. Equation 3.53 can be simplified further, to facilitate 


computer prcgramnung, 


ms 


be; SCCe, Se ee (3.54) 
; j=! 
and AT,, is defined as AT,, = CT, 
Al; =CT,,, + CTa,, J}=2.3.....N 


ATi v41 = Cie. 


In convenient matrix, the tangential velocitv can be written as 


Vi, cos( Fy a a} ATi) ATi 2 cee Any ATi v4) ms 
\ oe cas| As - a) ; AT; Atove, a 

= Oe . ; ( J.0on 
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5. Calculation of Cp.CL.CD and CNI 
After the velocities at control points have been determined. one can obtain the 
pressure coefficients at the i-ch control point using the Bernoulli equatien 3.37 
Similarly, other coeflicients can be obtained ‘castl by. pressure Intecrumeonma ome meme 


body contour. 


D. DISCUSSION OF THE PROGRAM PANEL 

The Panel program used in this analysis consists of two FORTRAN programs 
which involve the implementation of source and vortex distmbutions, and the 
implementation of linear vortex distributionsm Onis senemmist mrocramets imc Wdedeinecans 
thesis (see Appendix A) as 4 Sample procera meee im meme Ouse uscd ence in. 
basicaliv the same as the Hess and Smuth method [Ref. 4] where the source is located 
on the mid-point of each panel, consraint alone the since crentmEOmmpaiel to 
panel. On the other hand. the vortex is constiered@eonsi2 men one dome mecc!s sulla 
satisiv the Kutta condition at the trailingedge, the veloc acommencmumane commoucce 
by differentiating the velocity potential in the horizontal and vertical directton. And 
introducing these values into the Kutta condition, we can solve for the unknown vortex 


strength and suosequently for the source strength. 
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Source and Vortex Panel Method. 
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In the second method, only vortices are used as singularity distributions. The only 
difference here 1s, that the vortices are distributed linearly along the panel and 
continuous from one panel to the other. 
1. Input Data 
The input data for the program PANEL must be arranged in the following 
order: 
1. Header card. The header card consists of the input for the number of 
coordinate points (column 1-10, integer) and input for angle of attack (colunin 
11-20 real) 


Coordinates of the airfoil. The arrangement for the body coordinates must be 


ty 


inputed in the following sequence: start from the trailing edge, progress on the 
lower surface to the leading edge, return through the upper surface and finish at 
the trailing edge, so that the trailing edge coordinate will be accounted twice. 

2. Program output 


The output from this program (see Appendix B) can be arranged as follows: 


| _ TABLE! 
DESCRIPTION@P [HE PROG Kea eO et 


Te sieds: : is the panel number 
2 : x coordinate of control points (mid-points) 
Bel) : v coordinate of control points (mid-points) 
Ae Vite Teele) : is the dimensionless velocity for each control point* 


(total velocity divided by free stream velocity) 


5. GAMMA(I) : vortex strength(must be the same for the each panels) 
6. SIGMA(I) : source strength for each control point 

7 Ora ' pressure coefficients for each control point 

8. CL,.CD,CM _ : coefficient of lift, drag and moment respectively 


E. DISCUSSION OF THE RESULTS 
Figure 3.6 shows some results from the present computation for various angles of 
attack. In general. both methods give very good agreement in pressure distribution. 


although there are slight differences in the region of minimum pressure and at the 
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Figure 3.6 Comparison of Nees distributions on the NACA 23012 airfoil using 


the original Smutn-tless-pance! metnod «source and vortex panels) 
and the vortex punel metnod.. 
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trailing edge. In the leading edge neighborhood , the source strength may give a 
significant effect in the region where the velocity changes rapidly. On the other hand, 
the different implementation of the Kutta conditions at the trailing edge also give slight 
differences in the pressure distribution. Recall that the first method using both 
velocities at the first and last panels are the same in magnitude but in the opposite 
direction, while the second method is using the condition of zero vortex strength at the 
trailing edge. 

From the above conditions, the conclusions can be obtained from the present 
computations. Both methods, can be used to predict forces and moments around the 
airfoil as long as the fluid assumed remains inviscid and is not changing with time. 
Experience dictated that the second method gave less execution time than the first 
method might be true because in the second method solve onlv one unknown (jy) 
instead of two (¥ and o) for the first method. On the other hand, the second method is 
more complicated in the numerical ‘formulation than the first one, because of the 


varving vortex strength throughout the airfoil contour. 
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— IV. VISCOUS FLOW METHOD 


In fluids. momentum is transferred by internal stresses, namely the hydrostatic 
pressure and the viscous stresses. When fluids are affected only by pressure and not by 
viscous stresses, their behavior is relatively easy to predict by standard inviscid flow 
methods as described in chapter 3. 

A variation of velocity in the direction normal to the direction of the velocity 
itself is called a shear. Especially in high Reynolds Number flows this shear laver is 
very thin. 

The most common tvpe of a shear laver is the boundary layer between a stream and a 
solid surface. On the solid surface, the fluid velocity is reduced to zero (no slip 
condition), but there is no direct constraint on the velocity gradient at the surface. At 
the outer edge, the velocity tends aswmptotically to the free stream values. 

As mentioned before, in 1904 L. Prandtl came up with his well known theory of 
boundary lavers. Prandtl’s hypothesis divides the flowfield past a body into two 
separate regions, namely: 

1. The region very close to the body where viscous effects are important. 
ro) )Me Tremiainine region where imertia terms are more dominant than viscous 
terms, so that this region can be treated as inviscid flow. 
These assumptions allow us to deal with the parabolic boundary layer equations, 
instead of the elliptic Navier-Stokes equations. The prior experience indicates that 
parabolic equations can be solved very rapidly and efficiently. Numerically, the change 
of characteristics means a change from a field procedure to a marching method, which 
integrates the boundary laver equation for given initial conditions step by step, 
proceeding in the downstream direction. Depending on the boundarv conditions, 
boundary laver methods fall into three tvpes: 
1. The direct boundarv laver method. This method emplovs the ‘no slip 
condition’, requiring zero normal and zero tangential velocity at the surface, 


and a prescription of the external velocitv at the edge of the boundary laver. 


tJ 


The inverse boundary layer method. The prescription of wall shear or 


displacement thickness replaces the above edge boundary condition. 


3. The interactive boundary laver method. The edge boundary condition 
prescribes a combination of displacement thickness and external velocity. 
Further discussion Will focus on the first and the third method. since the computer 


code uses those. 
Ae DIRECT BOUNDARY LAYER METHOD 


Direct methods are used tn the region where the viscous effects remain small. The 
current code applies this method in the region near the leading edge. Direct methods 
allow the generation of initial conditions at the stagnation point and the efficient 
integration around the leading edge. The numerical approach features a finite 
difference method. which recasts the continuity and momentum equation as a svstem 
of linear algebraic equations [Ref. 13.] To begin with, wevcomsider o-by s¥cadviiows oF 
incomipressible fluids. described in a curvilinear coordinate svstemi with x directing 
along the airfoil surface and v perpendicular to the airfoil surface. The velocity 
components uw and ¥ shall be determined such that they satisfv anywhere in the flow 


field the continuity equation (4.1) and the momentum equation (4.2) 


Ou Ov 
segue i (4.1)- 
gg On| 2 as Oe a Oneal gc 
Oz dy =e ' dy Oy rs) 


where 3 boundary conditions are needed at the boundaries of the flowfteld, 
y=: i wise) ei). Be. 0) eau 


Y= Yet ult ye) = ue(2) 


With bi: clases vv. These equations are referred to as boundary layer or thin shear 
layer equations. To solve these equations, It 1s convenient to introduce a stream 
function (u = Gy éy and v = - Gw ox ), which reduces thesmmmieer cderentent 
variables. Since the stream function automatically satisfies the continuity equation, 


only the momentum equation 1s left 


Ou Ore Ow O%n di, 0 Oru B 
Pa ae ee | a (4.3) 
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This equation 1s subjected to the Falkner-Skan transformation. which scales the normal 


coordinate y and the stream function y with reference to the external velocity. 
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With these transformations, the momentum equation beconies 
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and the boundary conditions are 
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where » 1s defined as a dimensionless pressure gradient parameter { om 
and prime denotes differentiation with respect to yn. This is a third order partial 


differential equation. -The solutions of this PDE are called non-sinular flows, because 


’ 
5 | e 


rf 


(4.4) 
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thex depend on both x and n. In contrast. if the right hand-side of equation 4.4 


teeemes, anc therefore the solutions depend on y only, they are known as similar 


flows. 
1. The Box Method 


0 


One of the inmost flexible methods in solving non-linear differential equations 1s 
the box method developed by Aeller (Ref. 13} 1n 1970. The basic steps of the box 


fiewvoGeate the canversion of the governing equations into a first order svstem. the 


discretization of the differential equation by using central differences and two point 


averages, the linearization. and the solution of the resulting algebraic svstem. 


The 


introduction of two additional dependent variables U and FV converts the third order 


momentum equation (4-4) into a first order svstem 


fi=T 


(4.5) 


(4.6) 
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and the boundarv conditions become 


ys 0e Ce) =0) f(2.Q} = 0 

R=%e: Ule.yj=l 
Instead of dealing with continuous functions f, Ul, and J’, we use a set of discrete 
values of these flow properties. Let the soluucn domain Qi= 2 = 2a) eee 


covered bv a rectangular mesh (Figure 4.1) 
oO) fj=2,-,-% with 2<i<f 


ny = 0. f= ep ch, with See aed —ae 
We approximate all quantities whether it 1s a function or its derivative or a parameter 
like ‘b’, in terms of nodal values and coordinate of the network. The stream function 


and its first and second derivatives with respect to Ware aperevidied at Wie nmedecmes 


network by 
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The solution of the parabolic boundary laver equation at a certain streamline posiuon, 
sav x.. depends solely on the solution of upstream positions say x;7) .X;,>.-. while 
no downstream influence has to be considered. The overall solution can be obtained 
step by step with the-calculation propagating from the stagnation point into the 
downstream direction. The advantage of using first order equations and central 
differences is that we can reduce the domain of dependence from all upstream x- 
stations to the imimeédiate preceding one, hencevone stem Gi siiemsel tions proce dine 
writes the governing equations for a column of net rectangles (boxes) residing in the 


subdoniain 
we, S ee rc; ane 0 


and solves subsequently for the nodal values of the downstream face of the rectangular 
shaped subdomain. The x-station being currently solved holds therefore the superscript 
indicated by the term central differences the equations are satisfied midway between the 
nodes. 

The two ODE 's are centered about (4). 7 )_~4,) 


And the PDE is centered about (+,24/2-4)-1/2} 
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Figure 4.1 Net rectangle for finite difference approximation. 
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So equations 4.5 .4.6 and 4.7 can be written in terms of finite differences 
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and the boundary conditions take the form 


Ui =) fi =0 
a 
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The subdomain under consideration consists of J-/ net rectangles, the flow quantities 
of each being related bv a momentum equation. The equations 4.8 and 4.9 link the 
dependent variables to their n-derivatives. | 

-. Newton’s Nlethod. 

Unfortunately, the unknowns appear in nonlinear combinations. Therefore we 
introduce .Vewrton ’s Alethod to solve this nonlinear svsterm The solution involves an 
iterative procedure, in which the variables are linearized around their + aiues oueuie 
preceding iteration 
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where k denotes the iteration counter. Substitution of these values into eee 4S, 
4.9 and 4.10, and dropping the quadratic terms in ea wo ee) lead to a 


linear system in the unknowns eae va ee 
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where the coefficients in the bracket are defined as 
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3. Keller’s Block Elimination Method 
Together with the boundary conditions this svstem is repeatedly solved until 
6fe* su BV —are small enough to be neglected. Equation 4.11, 4.12, and 4.13 
can be solved by Acller’s Block Elimination Method. Biock-tridiagonal matrices are 
composed of submatrices. called blocks, of which only those residing on main and both 


adjacent diagonals have non-zero entries. 
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and where the righs hand sides are obtained Som 
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The unknowns of the linear equations are the Newton iterates of the stream function ( 
Ts ), 1ts first derivative ( Fie ) and its second derivative ( Be ). This method is 
very effective and it consists primarily of two steps: 
1. The forward step eliminates the lower diagonal of submatrices. 


2. The backward step solves the remaining system from bottom to top. 


B. INTERACTIVE BOUNDARY LAYER METHOD. 

The application of the direct method is restricted to regions where the viscous 
effects remain small. Integration of the boundary laver equations will break down at 
the point of zero skin friction. To avoid this break down, we need a method that is. 
able to integrate the boundary laver equations through the point of zero skin friction. 
PUeederetnese methods are required to account for strong interaction effects, arising 
from boundary laver separation or the rapid flow acceleration downstream of the 
trailing edge. both of which cause substantial changes in the external velocity 
distribution. : , es 

In contrast with direct and inverse methods, the interactive method treats the 
external velocity and displacement thickness as unknown quantities, reflecting the 
elliptic character of the outer flows. This introduces apparently one additional 
unknown into the viscous flow problem. whose solution can be obtained by using two 
methods : 

1. The eigenvalue method, or 

2. The Mechul function method. 
The second method is being preferred here, since the first method involves non-linear 
eigenvalue problems. The edge boundary conditon of the direct problem 1s 
supplemented bv the so-called interactive boundary condition, which relates the 
unknown external velocity with its inviscid and “displacement-perturbation-related 
contributions. Boundary laver equations in the following constitute a system in the 


unknown functions u{x.y), (xy) and u (xy 
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These equations consist of continuity equation (4.15) , momentum equation in x-direction 
(4.16) and the seemingly unnecessary omentum equation in y-direction (4.17). where 
the pressure term has been exoressed in terms olnicee semaine ecrociine The Mechul 
function approach assumes that the external velocity be a function of two argunients, 
resulting in the need for the trivial v-momentum equation. The reason for considering 
u(xy') rather than w(x) is for purely numerical reasons. 1e., such a provision allows an 
easv setup of the finite difference equations avoiding the eigenvalue technique. 

The govering equations are complemented by proper boundary conditions. The 
velocity components u and v are required to satisfv the no-slip conditions at the surface. 
of the airfoil and the horizontal component must merge smoothly into the outer flow 


at the boundarv laver edge. 
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with ux) denoting the inviscid velocity distribution and the second term, called Hilbert 
integral, approximating the perturbation velocity due to viscous effects. The interactive 
method can be applied to attached and separated flow regions, while direct method 
cannot do so because of their breakdown related to the Goldstein singularity, nor 
inverse methods because of their pOor cOnvercence Miatesemne ne lorest mem erict) c 
methods are being preferred on the main parts of the airfoil. Only at the stagnation 
point this method cannot be applied. 

The steps which turn the parual differential equations of the interactive problem 
into a linear svstem of algebraic equations resemble those of the direct method, so that 


only major steps and their results will be repeated here. -After the introduction of a 


Stream function (u=Cy Cy and v =-6y.dx), the equations undergo a transformation. 
Which scales the normal coordinate v, stream function y. and the external velocity u, 
with reference to the constant velocity uy and the local streamwise coordinate x 


tbs) 





=) ae 
Vr 
i | 
AS) Sa 
Jf lye 
- Uel 5. ¥) 
ropapappeaniieicalll 
Un 


Where up 1s taken as the free stream velocity. The concept of constant boundary layer 
thickness, attained by Falkner-Skan variables with u, as reference velocity, has to be 
abandoned because the external velocity 1s wiknown in the interactive calculations. 
Provided that the integration of the boundary laver equations does not start in the 
Memnicdiate neichoorhood of the stagnation poimt. the growth of the boundary laver 
thickness can ode kept linuted. In terms of these so called sesmu-transformed coordinates 
the boundanm laver equations. written as a first order svstem by means of two 


adJiuional dependent variables U and J’. and the boundary conditions take the form. 
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The discreuzation of the flow field follows closely the above outlined procedure of the 
direct method, covering the generation of an orthogonal grid and the introduction of 
central differences and two-point-averages. The overall solution proceeds in the 


downstream direction, accounting for downstream travelling disturbances only. On the 


assumption that backflow velocities are comparatively small, a stable integration can 
be carried out by adopting the FLARE approximation (Flugge-Lotz and Revhner). The 
purpose of a FLARE approximation ts to pernut the use of a downstream-marching 
algorithm in regions of backflow. This 1s accomplished by setting the streamwise 
convection term u Cu Gx equal (6 zero im fee-ons ec) @aanino ww, Wich 
1 ee Gt 
De yo : Umi? 29 
J~lYf2 '¢ [7s 
0 rae j-2 < 9 
designating an “on-off switch” of the streamwise convection term, the finite difference 


equations of the interactive boundary layer problem become 
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Boundary conditions are expressed in terms of nodal values. whereby the evaluation of 
the integral occuring in the interactive boundary condition involves an approximation 


in the fashion of the panel method approach. leading to 
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where z, and c.. denote a parameter and the diagonal element of the interaction matrix, 
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resulting from a discrete approximation to the Hilbert integral. Averaging as well as 
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centering 1s supposed to obev the principle that the number of generated terms 
approaches a minimum. This entails ordinary differential equations, like the y- 
momentum equation, being centered about the middle of the downstream face and 


partial differenual equations, like the x-momentum equation, being centered about the 


middle of the box. 


A balance of unknowns, which occur here as vectors in four components. 


fv), 171, }* confirms the principal solvability of the svwstem. The J quadruplets of 
unknowns match with a total of 4/ equations, including 2(J-/) auxiliary relations, (J-/) 


peat, 


X-momentum and (J-/) v-momentum equations, each of which corresponding to one of 
the (J-/) net rectangles. and 4 boundary conditions. After linearizing this svstem 
Gemmerthe values of the preceding iteration (iteration counter “K-/"), respectively 
arOund the solution of the adjacent upstream x-station in case of the first iteration, we 
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supplemented bv the two components of no-siip condition. edge and interactive 
oe@@ndar. condition 
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Terms have been grouped such that known quantities reside on the right hand sides, 
while unknown quantities appear left of the equal sign. The abbreviated coefficients in 


the momentum equation are defined by 
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and the right hand side of equation (4.28) 
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Since the overall procedure involves a repetitive linear pattern to approach the solution 
of the nonlinear svstem. the linear equation solver has to be as fast as possible. Fast 
algorithms take advantage of specific matrix structures. one of which is the block- 
tridiagonal structure. which can be enforced in this method bv the following 
arrangement of equations and boundary conditions : 
1. First set of equations (indexy = / ) 
|. Wall boundary condituon prescribing no penetration 


ape 


2. Wall boundary condition prescribing zero tangential velocity 


Gz 


3. Second auniliary relation (linking U” and V) of the first box 


+. Trivial y-momentum equation of the first box 


tJ 


Intermediate sets of equations (index 2 SjS/J-/) 

1. First auXiliarv relation (linking f and U) of the (j-/)-sz box 
2. Momentum equation of the (j-/)-sr box 

3. Second auxiliary relation (linking U’ and V) of the /-ch box 


~. Tnvial v-momentum equation of the j-rh box 
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Last set of equations (index j= J) 

1. First auxiliary relation (linking f and U) of the last box 

2. fomentum equation of the last box 

3. Interactive boundary condition 

“. Edge boundary condition. 

Following these instructions equations and boundary conditions can be arranged in 


miaterx-vectcr form. 
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(rayes}F are four dimensional vectors of the unknown Newton iterates and the 
known right hand sides. respectively. The blocks in the three diagonals of the above 


Matnx ate —x4 and can be obtained from 
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The components of the right hand side vectors can be calculated from the following 


formulae 
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The numerical solution of the above svstemimcam deaim cemaemle comeum ems 
olock elimination method. which works very much like Gauss's algorithm. but. firstly, 
matrices are eliminated instead of scalars, and, secondly; so Uiteeemie eine et omosean 


be saved because of sparse occupation. 
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C. INTERACTION MODEL 


The interaction model refers to the coupling of the boundary laver and the 
external inviscid flow. In principle, the interaction consists of thickening the effective 
airfoil shape by viscous displacement. which will result in a change of the surface 
pressure. Numerically the interaction between viscous and inviscid flow regions takes 
place via the boundary conditions only. which is the specification of either an 
impermeable displacement surface or a nonzero wall transpiration at the original 
surface in case of inviscid flow, respectively, the prescription of either pressure. or 
displacement thickness. or a linear combination of both in case of viscous flow. If the 
viscous effects on pressure remain small. then the interaction is called weak. However. 
situations, where viscous disturbances to the inviscid flow field are substantial, demand 
the application of strong interaction. In general , interaction models can be classified 
into four types. whereby the first three types provide loose coupling of viscous and 
inviscid regions: 

|. The direct interaction method combines a direct inviscid and direct viscous flow 
solver. This classical approach achieves a solution by iterative matching of 
boundary laver calculations with prescribed pressure and inviscid calculations 
with prescribed displacement thickness. The alternate treatment of pressure and 
displacement thickness leads to a local breakdown of the procedure when slight 
changes in pressure entail significant changes in displacement thickness (see 


Figure 4-2a). 


tv 


In the inverse method the roles of the displacement thickness and pressure are 
interchanged. Hence the inviscid flow equations deternune the displacement 
thickness distribution, which is imposed as boundary condition on the viscous 
flow calculation. The result of which is the pressure. which closes the cycle bv 
being input to the inviscid flow computation. The hierarchical manner of 
solving the complete flow problem exctudes this model just as the previous one 
from handling strongly interacting regions (see Figure 4-2b). 


The semi-inverse interaction method is composed of a direct inviscid and an 


Gs 


inverse visous flow solver. Both parts of the scheme process the same input, 1.e. 
displacement thickness, and generate the same output, i.e, pressure. After each 
cycle an updated displacement thickness is relaxed based on the deviation of the 


two pressure distributions, which should coincide upon convergence (see Figure 
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Figure 4.2 Direct, Inverse and Semi Inverse method. 
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4-2c). The existence of a definite hierarchy between boundary laver and the 
outer flow in the above mentioned methods gives rise to a limited rate of 
feedback between both regions. Strong interaction models incorporate the outer 
flow somehow in the boundary laver calculations, for example by the following 
interaction law: 

4. Simultaneous or strong interaction methods solve the boundarv laver equations 
Subject to an interaction law as outer boundary condition, which retrieves the 
elliptic character of the outer flow. No definite assignment of displacement 
thickness and pressure can be made to viscous and inviscid region. Rather. both 
quantities are treated as unknowns, related by the interaction law. The 
procedure emphasizes simultaneous solution for both displacement thickness 


and pressure (Figure 4.2d). 


id 
The current interaction law is formulated in terms of displacement thickness 0 and the 
external velocity w,. the latter being related to pressure dy Bernoulli's equation. The 
SOolucion of the boundary laver equations in an iterative fashion makes use of an outer 


boundary condition. in which the tota!] external velocity is written according to 


e(r) = Uez( 2) + Ues( 2) (4.31) 


Where wx) is the inviscid external velocity and wig (x) is the perturbation due to 
displacement effects. Using the thin airfoil concept, the perturbation velocity can be 


Written as the Hilbert integral 
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The contribution due to viscous effects can be evaluated bv means of the blowing 
velocity concept. Lighthill proved that the effect of boundary lavers on the outer flow 
can be represented by a surface distribution of sources. The source strengths must be 
deternuned such that the virtual displacement surface becomes a streaunuine (see Figure 
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VISCOUS FLOW SOLVER 
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Figure 4.2d Viscous-Inviscid Interaction Method. 
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This implies that the blowing velocity v(.x.0) equals half of the local source 
Strengue 
=~. The displacement thickness is assumed toy Gemso smidtl that tie momen 


velocity components of the inviscid flow do not varv across the boundary laver 
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with(¢,,,;denoting a matrix of interaction coefficients defining the relationship between 
the boundary laver thickness and the external velocity. Recalling that the boundarv 
laver calculation at a streamwise position involves Newton's iterations, the inviscid 
contrioution can be included in the total external velocity of the previous Newton's 


iteration, leading to a generalized version 6! the imlenacr@ asian 
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Since displacement thickness does not belong to the dependent variables, (6")'*, the 
displacement thickness of the streamwise position currently being solved, must be 


exoressed in terms of dependent variables to make allowance for its unknown status. 
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Figure 4.4 Application of the Direct and Interactive Method. 
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With (6*)''* being replaced by Yaesu jer and after separating known and 


unknown terms. the interaction law takes the form 
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Implementation of this relation necessitates a known right hand side, which can be 
evaluated only in approximate fashion, because the term {u,5")*** is not not known 
yet downstream of the current x-location. To ensure interaction also over this region, 
these terms are taken from the previous iteration updated by some relaxation formula. 
Sal gee J ve / tty and = gS = 3"/u denoting the dimensionless interaction 
coefficient and right hand side, respectively,.the actually coded interaction law can be 
Written in terms of semi-transformed coordinates | 


x 


a (1 — Ci; ny) + Ci i = dj (4.38) 


This equation is being used as an outer boundary condition in the viscous flow 
computation, and relates the unknown external velocitv with the unknown 
displacement thickness of the i-th boundarv laver station, whereas the displacement 
thickness has been expressed by the streamfunction and external velocitv. Because of 
the elliptic character of the outer flow, which has been incorporated in the boundary 
laver via the interaction law, the solution requires a global iteration, consisting of 
several viscous sweeps to be performed over both the streamwise upper and lower 
surface. 

Let the continuous function “external velocity times displacement thickness”, 
denoted by D, be discretized at a finite number of streamwise positions, then the thin 
airfoil integral can be approximated by a finite series of weighted Ds at the very 


locations 
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The weights are the interaction coefficients c., , with the first subscript indicating the 
streamwise position. where the correction term us is to be evaluated, and the second 
indicating the location, whose effect is accounted for. With a properly interpolated D- 
function, integration can be carried out analytically piece by piece. Provided the point 
under consideration does not fall within the limits of integration, D will be interpolated 
linearly. AA piecewise linear function can be built up bv overlapping triangular 


distributions, integration over which yields the coefficients whose kzi, kzi-] , ke it] 
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A linearly interpolated D-function would lead to singular integrals for the coefficients 
K=1,x=i-1 andk=i+/. Therefore D will be approximated by a polynomial of degree 
2 in the interval x, , S 3 = x;.). Splitting again into overlapping distributions, 
which this time are parabolic, and applving Cauchy's principal value technique permuts 
integration of singular integrands. The coefficient at the middle of the inducing source 


distribution 1s given bv 
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The coefficient at the right of the inducing source distribution is given by 
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The coefficient at the left of the inducing source distribution is given by 
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As indicated above, the overall solution 1s approached in an iterative process, in 


which alternately viscous and inviscid flow equations are being solved: 


iE 


tJ 


Calculate the external velocity distribution w,, in an inviscid flow held by means 
of the conformal mapping method. To account for the airfoil-thickening due to 
viscous displacement, specify a nonzero normial velocity (blowing velocity) at 
the airfoil surface. 

March through the boundarv lavers of both streamwise upper and lower surface 
using the interaction law as outer boundarv condition. 

Check convergence and quit if the criterion satisfied. 

[f the convergence criterion 1s not met, prepare for anothep cvcie. Update the 
product-term “external velocity times displacement thickness” on the basis of 


the deviation between inviscid and viscous external velocity distributions 
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where A denotes the-counter of global iterations. Further, compute the blowing velocity 
( J, = Wall transpiration velocity ) distribution, which serves as boundary condition 
for the inviscid flow solution 


(ee \! = os; pot ose (4.46) 


and proceed with the first step. 


D. TGRBULENCE NIODELLING 

The precence of additional unknown shear stresses in turbulent flows requires 
modelling assumptions to balance the number of unknowns and equations. Eddv 
EeOSieh: models. one of which is used in the mrescnemmethod. [Ref. 13.) relate turoulent 
shear stresses to mean flow quantities on an empirical basis. Thev draw their versatility 
moo (ie-convenience of maintaining the same approach and numerical formulation for 
Det: lanunar and turdulent flows. sAccording to this formulation for wall boundary 
laver tlows. the turbulent kinematic eddy viscosity is defined bv two separate formulas. 
@g-8 8, the inner rézion being vased on Van Driest’s approach, and the other for the 


outer region being based on a velocity defect approach 
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Conunuity of the turbulent kinematic viscosity 1s established by defining 3. as the 
distance from the wall where expressions for inner and outer region do agree. Since 
transition 1S not an instantaneous process. an intermediate status of flow 1s assumed 
between the laminar and fully turbulent regions. This region is taken tnto account by 
Intrcducing an internuttency factor, which smears out the step-shaped change from 


kinematic to eddy viscosity. The formula here 1s suggested by Chen and Thyson. 
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where R,,, denotes the transitional Reynolds Number ( u,x/v),, ie., Reynolds 
nf ? 


Number based on external velocity and streamwise location x,, at onset of transition. 
A value of 1200 was originally assigned to the empirical constant Gorey: Flower ere 
numerical experiments seem to indicate that low Reynolds number flows can be better 
modelled by values below 1200 (further discussion in the next section). The parameter 


ad in the outer region formula is obtained from 


a= = - = 0.0168/[1- 9 Saal 


where the nondimensional factor F denotes the ratio of total turbulence energy 


production to the shear-stress-related turbulence energy production, evaluated at the 


location of maximum shear stress 
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The ratio of the time-averaged quantities above can be approximated by a function of 


Rr => ty /( we) maz 
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Since this turbulence model is not validated for separated flow, eddy viscosities in 
regions of backflow correspond to those of adjacent upstream station at the same n- 
coordinate. Transition, if not available from other sources, currently is predicted by an 
empirical data correlation expressed in terms of the Reynolds Numbers based on 


momentum thickness and streamwise coordinate at onset of transition 
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E. DISCUSSION OF THE COMPUTER PROGRAM 
1. Inviscid Flow Nlethod 

The inviscid flow method used in the interaction method is based on conformal 
mapping Which can be divided into the transformation of the region outside the airfoil 
Semper reoion Oltside a wumit circle amd to the solution of the equations in the 
transformed plane. The transformation was achieved in three parts which together 
reoresent the major computational effort. 

I the first mapping. the airfoil is perturbed slightly to make the upper and 
lower surface trailing-edge points coincide. This is accomplished using a logarithnuc 
mi@epine function and is necessary onlv in those cases in which the airfoil trating edge 
has non-zero tnickness. 

In the second mapping. the trailing-edge corner is analvtically removed bv 
apolving tne Karman-Trefftz mapping. 

In the last mapping, the resulting quasi-circular shape is mapped to a perfect 
Girele wSing an iteruted sequence of applications of the fast Fourier transform 
alcorithm. The calculation of the flow in the transformed plane also makes use of 
meucereanalysis technique. . 

The major computational effort required in the inviscid flow method is due to 
Ure iransiemmation. lwisenecessary to compute thestransformation»soniveonce in the 
Viscousanviscid flow interactions. so that the computational expense due to inviscid 
calculations can be held to a minintum. When the angle of attack increases, care must 
be taken since the displacement thickness can become fairlv large, approaching 10°% of 
the airfoil chord. In this case, use of the blowing velocity ( equation 4.33 ) on the 
airfoil surface produces a dividing streanuine from the leading-edge stagnation point 
which approximates the edge of the boundarv-laver displacement thickness. The 
inviscid flow outside this dividing streamline is therefore the same as that pasi the solid 
body defined by this streanuine. However. inside this dividing streamline, the inviscid 
flow is fictitious. In particular, the assumpuon that there 1s no pressure variation 
across this ficitious region becomes invalid as tne magnitude of blowing velocity V’..,, 
imereases. [he approach adopted here is. therefore, to evaluate the velocity distribution 


directly on the displacement surface while still applving the blowing velocity on the 


original airfoil surface. To avoid a discontinuity of the velocity between upper and 
lower surface, the Kutta condition is applied on the dispiacement surface. Requiring 
equal off-bodv pressures for the upper and lower trailing edge points, a quadratic 
equation can be solved for the unknown circulation. 
2. Interactive Viscous Flow Method 

In this method. the boundarv laver is transformed into Falkner-Skan and 
subsequently senu-transformed coordinates. The numerical solution of the transformed 
equations is performed using an implicit finite-difference technique, which 1s first-order 
accurate in the streaniwise direction and second order accurate in the normal direction. 
In principle, the computer program calculates: 

1. The inviscid pressure distribution for a specified angle of attack 

2. The boundary laver, including the velocitv profile, skin friction, momentum and 
displacement thickness either: 

1. subject to a prescribed pressure distribution, or 

2. subject to an interactive edge boundarv condition. In this case. the external 

velocity will be a part of the boundary laver results. 

These calculations are performed for a specified number of x-stations on the airfou and 
in the tvake, a number of sweeps is made on the airfoil in order to obtain a converged 
solution. The Cebeci-Snuth two laver model is used here for computing the eddv 
viscosity Vv, and the transitional flow region is modelled using- equation 4.48. The 
method also employs a semi-empirical formula to predict the onset of transition. This 
criterion, equation 4.51. was proposed by Michel. The code offers two possibilities how 
to deal with transition: 

1. The loci of transition are calculated bv the code. In this case Michel's criterion 
is eniploved to predict the onset of transition. When laminar separation occurs 
upstream of the calculated point of transition, then Michel's Criterion 1s 
disregarded and the onset of transition is redefined at the point of lamunar 


separation. 


tv 


The begin of transition is specified bv the user (fixed). The code has been 
modified to allow the onset of transition to be within or downstream of the 
laminar separation bubble. The previous version of the code always redefined 


transition at the point of lanunar separation. 
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F. DISCUSSION OF THE RESULTS 

1. High Reynolds Number Flows 
The computer program was first applied to high Revnolds-Number flows about the 
Wortmann airfoil FX 60-126. Figure 4.5 shows results for Revnolds Numbers from 
700,000 up to 2000,000. It is seen that the lift predictions are in quite good agreement 
with experimental data. 

2. Low Reynolds Number Flows 

The subject of low Reynolds Number flows is important to a number of 
practical problems, including remotely piloted vehicles, wind turbines and propellers, 
sail planes. human powered vehicles, etc. Viany boundary laver phenomena which have 
eluded accurate analvtical prediction occur within this flow regime and are associated 
with laminar separation and subsequent transition of a laminar free shear layer. Figure 
4.6 shows the phenomenological features of the boundary laver on a low Reynolds 
number airfoil. a 
When the Reynolds Number based on momentum thickness, Rg, is sufficiently low, the 
boundary layer remains laminar up to, including, the point of minimum pressure or 
maximum suction. At some location between the nunimum pressure and the 
theoretical point of laminar separation the -Revnolds Number of the boundary layer 
attains a critical value. This is referred to as the transition Reynolds Number based on 
momentum thickness, Rorr | 

The provision of reliable information of the flow characteristics of low 
Reynolds Number airfoils'is hampered bv the sensitivity of the flows, and particularly 
those involving separation and transition. 

In recent vears, several theoretical studies have led to the development of 
methods for predicting airfoil characteristics at low Reynolds Numbers, but there are 
no universal methods which can accurately predict and account for a separation bubble 
in the design of efficient low Revnolds Number airfoils. 

The viscous/inviscid interaction method was applied to the NACA 65-213 
airfoit at a Reynolds number of 240,000. It was found that the calculation would fail 
to converge if transition was predicted by Michel’s criterion (equation 4.51) and if the 
empirical constant Guty was chosen to be 1200. Therefore it was decided to investigate 
the influence of the start of the transition and of the constant Garey on the resultssbm 
systematically varving both parameters. Table | of Appendix C shows the predicted 


lengths of the separation bubble which is obtained if transition 1s chosen to start at 64, 
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Figure 4.6 Phenomenological features of the boundary laver 
on the low Réevnolds Number airfoil. 
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Figure 4.7 Experimental result for the NACA 65-213 (bv Hoheisel et al.) 
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(a). Velocity distribuuon and (6). Flow visualizauon. 


Gomes) G7. GS, 69, 70, 72. 74, or 76 % of chord and if Gyty is chosen as 10, 20, 40, SO, 
or 120. It is seen that there are significant changes in the length of the separation 
bubble depending on the chosen parameter combination. This effect is displaved more 
clearly in Figure 4.11 and 4.12. 

An increase in Guty increases the transition length as well as the length of the 
separation bubble. 

Hohetsel et al. ([Ref. I4] performed detailed laser- Doppler velocimetry 
measurements of this airfoil at a Revnolds number of 240,000. Their results are shown 
fo tieures +7, +9) and 4.10 and at was attempted to chooSe that parameter 
combination which would give the best agreement with the experimental results. If the 


begin of transition ts chosen at 74 °% of chord and if G.,,.= 20, then the results shown 


a AOg 
in Figures 4.13, 4.14 and 4.15 are obtained. It is seen that the boundarv laver profiles. 
displacement and momentum thickness distributions upstream of the separation bubble 
are in excellent agreement, Whereas considerable deviations are found in the bubble. 
Finally, Figure 4.16 through 4.20 show the computed boundary laver velocity profiles 


for different values of Gutr ranging from 10 to 30 
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Figure 4.9 Effect of variation of NTRUL for Gay = 20 


on the velocity distribution. 


U/UE 


U/UE 


GGTR=40 GGTR=40 





CGTR=40 
XTRU=0.74 





02 


_ XTRU=0.69 XTRU=0.70 





CCTR=40 
XTRU=0.76 
0.4 0.6 0.8 1.0 2 0.4 0.6 0.8 1.0 
PERCENT CHORD(X/C) PERCENT CHORD(X/C) 


Figure 4.10 Effect of variation of XTRU for G. 
on the velocity distribution. yes 


86 


U/UE 


GGTR=40 
ATRU=0.72 


~ 





A) a2 0.4 0.8 0.8 
PERCENT CHORDXX/C} 


VELOCITY DISTRIBUTION 


‘NACA 65-213 


REYNOLDS NO:240,000 
ANGLES OF ATTACK:0 DEG 


NOTES 


GGTR: EMPIRICAL CONSTANT 
XTRU: BEGIN OF TRANSITMON(UPPER SURFACE) 


40 


Skin Friction Coefficient, C, *103 


Intermittency Factor, Pay, 


0.0 


| 
O 


. 
ae 
O 


-3.0 


ar 
Bie S| a] ks ia VI ie | ty T | biel per [ Vue | a oe | 


a 
O 












AIRFOIL DATA 





fF NACA 65-233, 
L“- Angie Of Attack: O ceg, 
~ = Reynolds No.: 240000, 
™. Start Transition: 0.68 
a 
: G 
On SN G_,.=20 
E | (Ax/c),»= 0.12! i. 
; WI tae eo is Gy, =40 
L its. (xver-— ©.15 
- $$ $$$ $e G,,,=60 
| (Ax/c),,= 0.18 
: G_,,=80 


! { ' ete eatin t 


00.55 0.60 0.85 0.70 0.75 0.80 0.85 0.90 0.95 100 
Chordwise Distance, x/c 





OC 
un 







0.09 
(Ax/c)eg= 0.14 
Ax/c)+,= 0.20 


| ‘ | Ax/C)en= 0.26 | 
| 


— 


O 
Qn 


Cherawisce Miscance, x/c 


cae 4 ll Effect of vabianion of G... for \TRL =1.68 


on tne distributions of skin friction and Untermuttency factor. 


87 


ero oeTOe morc (om o womee oO. 0.85 0.90 0.95 1.00 


5 
Oo 
















foe) = 
Se AIRFOIL DATA 

a NACA 65-213, 

C) 2.0 Angle Of Attack: O deg, a 

= Reynolds No.: 240000 

S 10 Guy au 

= 0.0 Stee =F 

© | RY 

CUE WA) | : 

SC -2.0 ee IN (x/c),,=0.66 
‘S (Ax/c),,= 0.12 Ta oc 
= “eee (Ax/c),,= 0.16 | (x/e)..=0.70 
is ! n=O. 
Cc -4.0 ___'Ax/c)s,= 0.21 (x/c),,=0.72 
Fe BS 0) (Ax/c)..= O25 (x/c),,=0.74 

-6.0 : 


0.500.55 0.60 0.65 0.70 0.75 0.80 0.85 0.80 0.95 1.00 
Chordwise Distance, x/c 






(Ax/c}ra= 0.14 | 
(Ax/c)rya= 0.14 

Ax/Clra= 0.15 

L J (&x/chin= 0.16 








(x/c),,=0.74 
| Kf (x/c},,=0.72 
| . (x/c)4,=0.70 
vat 








Intermittency Factor, 7, 
© 
'@ 4) 


(x/c),,=0.68 







oe a 
e) — 
“US a i i Sn i 





0.500.55 0.60 0.65 070 075 0.80 085 090 095 100 
Chordwise Distance, x/c 


Figure 4.12. Effect of variation of NTRU (begin of transition) for G. ie 
7 on the distributions of skin friction and internuttency factor. ! 


88 


3 





1 
O6 
1 NOSE 5-21 9 
D Re. 240.000 
N AOA=0 deg 
— oO 
S 
O38 
~— | GETR=20 
o | XTRU=0.74 
i> 
x x/c=0.367 
— 
© 
Co) 
See G.. 65-210 
O Re.240.000 
AO A=( deg 
<a) 
ce 
o 
ore 
> 
ml Cenk — 20 e 
~ NTRU=t).7 a 
= aa waa ——_ EXPERIMENTAL 
(e) 
S ss229+* INTERACTION MTD 
= 
nN ‘ 
o 
=| Nees os 
Re.240.000 
00 eae 
O AOA=0 deg 
S 
< GGTR=20 
oe XTRU=0.74 
© 
© 
= 
S 
O 





} } j 


00 02 04 06 O08 #10 
U/UE 


Figure 4.13 Velocity, profiles in front of the bubble 
aGk C = 07, US ane u,009. 


89 


v/C 


0.0000.004 0.008 


Whe 
0.008 


0.016 0.020 


OPO 


0.016 0.020 


OO 


0.000 0.004 


NACA 65-213 


REYNOLDS NO:240,000 
AQA = 0 DEG 


GGTR=20 
XTRU=0.74 





=0.15 0.05 0.25 0.45 ies 0.85 


NACA 65-213 
REYNOLDS NO:240,000 
AOA = 0 DEG 
CGTR=20. 

XTRU=0.74 





More 0.05 0.25 a. 25 6016S 0.85 


U/UE 
4.14 Velocity profiles in the bubble region 
BIBLE ot ce = OL uIRCIOWmn 


DSTAR ; THETA 


0 004 0.006 0.008 


0.002 


NACA 65-213 


“ 
me 
= 
REYNOLOS NO. 240,000 iS) 
ANGLE OF ATTACK O OEG rea 
EMPIRICAL CONSTANT 20 
BEGIN OF TRANSITION 4.74 = 
ay 
oy 
Y 
as 
re) 
wn 
Q 


INTERACTION METHOO --—-~- 


EXPERIMENT ----4 EXPERIMENT @-<-< 


INTERACTION METHOD -<< 


OZ 0.4 0.6 ORS 
PERCENT CHORD(X/C) 


Figure 4.15 Comparison of 6 and @ with experimental results. 


a 


el lee. 


seges — AIRFOIL DATA 


yic NACA 85-213, 
Angle Of Attack: 0.00 
0.02 Reynolds No.: 240000 
GGTR : 10 
SS 
0.01 -————__—_—_-_ — 
———— 





0.0 0.6 1.0 u/U,, —— 
Se oo! —— —_—— 
———_——. = 
Sey — > aay, 





Symbols: TI... Begin of Transition a ate 
T2... End of Transition ———— <a = 
ca —ae} 
SL.. Laminar Soparation ———— = 
ST... Turbutent Separation 4 —— = 
RA... Reattachment ———= — 
4.0 He IS 
3 
oe) 0.8 — 
3.0 oO 
= 2 
eal 0.6 
Oo 2.0 S 
5 ] 0 | — 
rs) Ss +. re : ; é 0.2 = 
= 3 
© ' = 
5 «8.9 6:0°-= 
O 
= ‘ 
Ono 7 
Ss ft 
Ue -2.0 - 
= \ 
ae \ 
YY 


-3.0 Vi, 


-4.0 





0.60 0.65 060 9.65 070 076 O80 086 090 0965 1.00 
Chordwise Distance, x/c 


Figure 4.16 Boundary laver profiles on the NACA 65-213 at Re = 240,000 
AOA = 0 deg and Gute = LV. 


on 


_ l,l FP + 


Figure 4.17 Boundary la 


Scales ~ AIRFOIL DATA 


y/e NACA 65-213, 
Angie Of Attack: 0.00 
0.02 Reynolds No.: 2400C0 
= GG7R - 20 


————_—_—~_ 

——— => 
— 
—— ——— 
=- 

0.00 





oe a = — - — - 
ee ee —— > ——s 
——————— Rare ——a ee eee =a —_——_—_ _ —— 
rt ——- a 
— 





Symools: Ti... Begin of Transition 
F2... End of Transition 
SL.. Laminar Separation ————— 
ST... TurDulent Separation 
RA... Reattachment 





4.00 - 
3.0 [ ; = 
ae —_— 
‘s ——e eee —— a! 
SS i (ee 
2.0 SS SS 
= / SS EE SS 
———— ee — ——! ———) ee ey 
0 SSS SS SS SS ESS 
- =ite ~ aa aa eee, lle ——ay . ead 
eee —_ — 
Sr ee Pie ( ae ee es 


Skin Friction Coefficient, C, *10° 
2 ! : 
Oo 
| 
i 


o) 
| (Le eu | (ate (i 
2) 
aes 
a 
D 
> 
—| 
NO 


0.60 065 0.60 0.65 070 O76 0.80 085 090 O86 1.00 
Cnhordwise Distance, x/c 


pooltlessonstihcs\ .GzG5-215 at Re = 


deo undeG.,,, = 20. 


Ver 
AOA ie az 


: 


— emmy 
00 05 10U/U, = 
: —————— ey 
——, ~~, 


eid 


Intermittency Factor, 7u 


0.0 


240.000 


Scales AIRFO!L DATA 


NACA 865-213, 

Angie Of Attack: 0.00 
Reynolds No.: 240000 
GGTR : 30 





————— — 
0.0 06 1.0 U/Us === 





aa 

-_— 

_—— = = 
Symboats: TI... Begin of Transition ee — 
T2... End of Transition ——— == 

SL.. Laminar Separation ee 

ST... Turbulent Seperation ey 

RA... Reattachment 5 oo 

eee eS 





SOS 
=e 
ee el = ——i 
Zoo ha ae atl 
————i . 
Oe Oo | 
eS eS a Ss See 
a a oe j 
1.0 Fs 








—— fl em onal — 
a et So 4 Ve eae 
Sg ey lg 
i (ee at ai a _ —_ t 
—— 


¢ 
— cll + eer , j 7. a , 
Eee = 
le 7 ——— ee 4 yd 





ea es 


Skin Friction Coefficient, C, *10° 


0.50 0.55 0.60 066 0.70 075 0.80 / 0.865 090 0.95 
Chordwise Distance, x/c 


Figure 4.18 Boundary laver profiles on the NACA 65-213 at Re 
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Figure 4.20 Boundary layer profiles on the NACA 65-213 at Re = 240,000 
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V. CONCLUSION AND RECOMMENDATION 


Cebeci's viscous inviscid interaction prograni was applied to the analysis of 
steady two dimensional incompressible flow past a NACA 65-213 airfoil at zero angle 
of attack at a Reynolds number of 240.000. Predicted boundary laver characteristics 
were found to be quite sensitive to the choice of boundary laver transition begin and 
length. Good agreement with the experimental results of Hoheisel et al could be 
obtained by proper choice of both transition begin and length. Further detailed 
measurements and calculations for other airfoils at low Revnolds number are 
reconimended in order to further validate the predictive capability of the 


VISCOUS. InViscid interaction method. 


a 


ae APPENDIX A 
FORTRAN PROGRANI 


CSNOEXT 


THIS PROGRAM CALCULATE THE PRESSURE DIS TRESUEO See eae 
AIRFOIL AT ANY ANGLE OF ATTACK, BY USING PANEL METHOD IN 
2-D INVESCrD, STBADYSELOn. 

IMPLEMBNPATION OF VORTEX AND SOUCRE DISTRIBU aie we eeu. 
ALL VELOCITIES AND LENGTH ARE NORMALIZED BY FREE STREAM 
VELOCITY AND CHORD™ LENGTH RESPe Gy Bile 


WRITTEN BY = CAPT.INDAP PHUTUT SUBROTe 
NAVAL POSTGRADUATE SCHOOL 
MONTEREY ,CA., OCIOBER ics 


KRREAKKKAKRKAKKRKKRKRARRRRRKRKRKAERRKRRKRKRKRKRRRKRKRKRRRRRRRRRRRRRRRRRRRERRKER 


OO) #4 AAA HARA HEED 
DS a, a SM Sa a, SS 


COMMON /NODE/ N 
COMMON /AAA/ AK(100, 100 ARK (100,100 ae eee (109-00 
COMMON /BBB/ BB(100,100),BBX(100,100),BBY(100,100) ,BBT(100,100 
COMMON /STAR/ BSTAR(100) ,CSTAR(100) 

COMMON /ALBE/ ALPA(100), BETA(100) 

COMMON /VEL/ VI(100) , VXI(100) , V¥I(100) 
COMMON /X¥/ X(100),¥(100) ,XB(100) , YB(100) 
COMMON /SSS/ $(100),SIGMA(100) ,SUMB(100) 
COMMON /ANGLE/ AN,ANG, TPI,TH(100) 

COMMON /PRESS/ CP(100$ 


C 
C ==--INPUL DATA 2 er NODE a Ones aa 
C 


READ(5,1) NN,AN 
1 FORMAT (2107 F10n) 


A} 


PI = 4.*ATAN(1.0) 
ANG = AN/180.*PI 
Venn = BS 
VYA = -SIN(ANG 
Rol = sore 

DO 10 I 


= 1,NN 
READ( 5S, 5) er 18 
10 CONTINUE 
5 FORMAT(2a1025) 


S 
RAERKKAKRKKAKKRRKRKKKKKKKKKRKKRERKKKRKRKKKKRKKKKRRAKKRKRKKKAKKRKRAKK 


x * 
x COMPUTE THE ANGLES AND LENGHTS OF EAGH@FANELS, * 
2 CALCULATE THE COORDINATE OF GON@ROL FPOllies ‘i 


KREAREKRKKRKRKERAKEKARERKKRREKRRRRRKRKAKRKEKRKERERRRRRRRRRRRAKRKRARAKKRE 


N = NN-l 
SS = 0.0 
Deut 1 = 1 
R(T) = S (xB t) + XB (T+1) 
y(I = ,5*(YB(I) + YB(I+tl 
TH(I) = ATAN2 (YB(I+1)-YB(I XB(I+1)-XB(I)) 
S (1) = : RT ( (HB(T+1) -XB(E) kKO + (YB(I+1)-YB(1I))**2) 
= + 


15 CONTINUE 


REE KREKAKAKKKKKKKKRKKKKKAKAKRKKAKRAERKKAKARKRKKKKKKKKAKKKRKRRRKERK 
x * 


- COMPUTE TIME ENDEPENDENT INPLUENGE (GenR Rema iais ~s 
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} 


23 FOR NORMAL, TANGENTIAL ,X AND Y DIRECTION s 
Pas 
RAAKAAKAKAKRAKKAAKKKKRAKARKKKKKARKAKKKAKKKAKAAKAKKAAKKKAKA KAR 


DO 25 I= 1.N 


mM 
C8] 
a9) 


Onto OwY 
1 a 


BBT 


EMDIF 
25 CONMENUE 


Co che pot ead Se eee | 


= er 4 
“~~ © gee OG 


C4, 


HH HHHe ~ 
CO QQ GQ Q~n, 


Se oe GaUauuuu, 


~ ~ = = = ~ 


wt 


a i gt es et” 


* 


THEN 
= Ono 
=) (0) 
= -0.5*SIN(TH(I 
= 0. 5“€OS Clint 1 
= 0Q.5*COS(TH(I 
= 0.5*Seign 
= 0.0 
ee OS 
oe - Reptegs op Satie) 
Uae (eet J) eZ 
HAR) 
sled) 
PL Sta (¥(I)-¥B(J))*COS(TH(J)) 
+ 15 ( J) “(Si 42. 8) 
SJ) =, BtAtS 
= TPI*(.5*C*F - D*G 
= TPi*(.5*Oaiees C*G 
= el SO OE \ aE + oe eee 
= TP is =e) AF - COS(TH(J))7G 
= ee 
= -AAX(I,J 
= oy 
= AA(I,J 


n 


* SETUP MATRICES ,SET GAMA=1.0, SOLVE THE SYSTEM * 
* BY USING GAUSSIAN. ELIMINATION WITH PARTIAL PIVOTING * 
n FOR TWO RIGHT HAND SIDES. s 
KIKI III KK IKK KKK KKK KKK RK KK RK KKK KKK KKK 
Gow : 
DOmeeer =1 -N 
SUMB(I) = 0.0 
Demoee J =1.N 
680 SUMB(I) = SUMB(I) + BB(I,J) 
eg = -SUMB(I) 
CSTAR(I) = -VXA * SIN(TH(I)) + VYA * COS(TH(I)) 
ie ay = BSTAR(I 
ABN I,N+2) = CSTAR(I 
68 CONTINUE 
CALL GAUSS(2) 
DO 500 I =1,N 
ALPA i = saan 
BETA(I) = AA(I,N+2 
oo CONTINUE 
C ---- MEET KUTTA CONDITIONS ,SOLVE FOR VORTEX STREGNTH(GAMA) 
: USING QUADRATIC EQUATION ---- 
le= 020 
ati = 0.0 
pe = 0.0 
fey = 0.0 
BX1 = 0.0 
=nO20 


BAN 
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Bay = - 0 


BYN 0 
DO 51 J =1,N 
AX1 = AX1 + AAX(1,J)*ALPA(J)+BBX(1,J 
AXN = AXN + AAX(N,J)*ALPA(J)+BBX(N,J 
AY1 = AY] + AAY(1,J)*ALPA(J)+BBY(1,J 
AYN = AYN + AAY(N,J)*ALPA(J)+BBY(N,J 
BX1 = BX1 + AAX(1,J)*BETA(J 
BXN = BXN + AAX(N,J)*BETA(J 
BY] = BY1 + AAY(1,J)*BETA( 
BYN = BYN + AAY(N,J)*BETA(J 


ane CONE ive 
EE = AXI“*~258¥ 0622 = elo eee 


PP AXI*BK1+AY1*BY1L - AXN*BXN-AYN*BYN +(AXN-AX1)*VXA 
+ + ({AYN@Avi)*VYA 
QO = BK1L**2+BY1**2-3XN**2-BYN**2 +2.*(BKN-BKL)*VKXA +2.% 
; 3 (BYN-BY1)*VYA 
R = PP*PP - EE*QO 
GAMA1 = (-PP + Se ae 
GAMA2 = (-PP - SORT(R))/(EE 
IF(ABS(GAMA1) .GT. ABS(GAMA2)) THEN 
GAMA = GAMA2 
ELSE 
GAMA = GAMA] 
: ENDIF 
: ---- SOLVE THE SOURCE STRENGTH SIGMA(J) ---- 
DO. S5e J =i 
Sieh GAMA*ALPA(J) + BETA(J) 
550 - CONTINUE 
Ee ; 
CALL CPRESS(GAMA) 
C . 
c . 
C ---- PRINT LIFT COEFFICIENT AND MOMENT COEFFICIENT ABOUT 
: LEADING EDGE ---- 
CEX = O80 
Chey == 20m9 
: CH 486= 02g 
DO” 110° ewe, 
a Sy = ee 
DY = YB(I+1) - YB(I 
CEX= CFXK + cP tt RADY 
FY= CFY - CP(I)*DX 
M = CM + CP(L)*(DX*K(1)+DY*Y(I)) 
110 CONTINUE 
‘ea =e ae ANG) - See 
e)9) = CFKX*COS(ANG) - CFY*SIN(ANG) 
WRITE(8,120) CL 
WRITE(8,125) CD 
WRITE(8,130) CM 
120° FORMAT(// felon CE = Fie. > 
125 EFORMATG/, s<, "GD =" aloes 
130 FORMAT(/,15X, 'CHLE =e Ss 


10. 
Peete "COMPUTATION COMP CETED! 
STOP 


END 
KARAKAK AR RKRKKKKK KKK KKK KKK KARA RRR RRR KKK KKK KKK KR KR 


SUBROUTINE GAUSS 
SOLVE SIMULTANEQUS EQUATION WITH TWO RIGHT HAND SIDES 
Bi GAUSS ELIMINATION WITR PaARGIAL FlveiiiGe 
SOLUTEONS STORES IN COLUMNS NEOUS I SAlIpRiEe S425e> 
MATRIX AA(NEQHUS ,NEQNS+2). 

(a ="COEFFICIENT OF AUGHENTED MATRIX 


BE oe a Se, ae. 
tA ADA A SE 


100 


“3 NEQNS = NUMBER OF EQUATIONS * 
i" ERHS = NUMBER OF RIGHT HAND SIDES ‘ 


We ek eee OL Rs Ou TOO) e.) Oa Ole @ ele ee 18. 16. ee 


SUBROUTINE GAUSS (NRHS) 


COMMON /NODE/ NEQNS 
COMMON /AAA/ AA(100,100),AAX(100,100) ,AAY(100,100) ,AAT(100,100) 
NP = NEQNS+1 


NEONS+NRHS 
GAUSS REDUCTION 
DO 40 I = 2,NEQNS 


SEARCH FOR LARGEST ENTRY IN (I-1)TH COLUMN 
ON OX BELOW MAIN DIAGONAL 


e 


NTOT 


OO” 


ANNAN 


IMAX = IM 
= ABS(AA(IM,IM)) 

D0 10 J = I,NEONS 

BMAX .Ge. ABS(AA(J,IM))) GO TO 10 


ABO( LE (J , mt) ) 


was 600 sq 
10 CONTINUE 
SeaeCH (i-1)TH AND IMAATH EQUATIONS 


AAAN 


Cy 
] 


Pe ISX .NE. IM) GO TO 30 
DO 


= IN,NTOT 


AA (Z 
29 CONTINUE 


ELIMINATE (I-1)TH UNKNOWN FROM ITH THRU 
(NEQNS)TH EQUATIONS 
30 DO 40 J = I,NEQNS 
R = AA(J,IM)/AA(IM,IM) 
DO 40 Ko aon 
AA(J,X)- = AA(J,X)-R*AA(IM,R) 
40 CONTINUE 


BECK SUBSTITUTION 


Bo. 70 a =98P ,NEOT 
AA(NEONS,K) = AA(NEONS,K)/AA(NEONS ,NEOQNS) 
DO 360 L 5 


AANA 


AAD 


= 2,NEON 
I = NEQNSt1-L 
I> = 
0 50 J =IP,NEONS 
50 Be a = ene eae 
60 REX IK =e x) /RA(I I 
70 CONTINUE 
RETURN 
=rod 


x 
* SUBROULSHE CPRESS 

a Se Ceesepeee oe essURE COEFFICIENTS AND TOTAL VELOCITY 
x 

x 


> a 


Ae TIED PORTS (COMPROL POINTS) 


Cr Sy ee ee i Oe es ee eee 6 6 6 8 cel we Og 8 68 te ee ee ee ee Ue ep ee eee 


SUS.OUTTNE CPRESS (GAMA) 


1OL 


exe) 


COMMON /AAA/ 33(185- 10) ‘BREi00; 100) BBY{ 200; 100), T(100, 100 
COMMON /BBB/ BB(100,100),BSxX(100,100),BBY(100,100),BST(100,100 
COMMON /PRESS/ CF (100 ) 

COMMON /NODE/ N 

COMMON /ANGLE/ AN,ANG,TPI,TH(100) 

COMMON /SSS/ $(100),SIGMA(100) ,SUMB(100) 

COMMON /VEL/ VI(100),VXI(100) ,VYI(100) 

COMMON /XY/ X(100),Y(100) ,XB(100) , ¥B(100) 


ao 


ae 
5 FORMAT (//, 15K, NEGA 23012") 
WRITE(8,85) N 
WRITE(8,7) AN 
6 FORMAT(15X,'NUMBER OF PANELS =',1I3) 
7 FORMAT(15X,'ANGLE OF ATTACK =!,F7.4) 
WRITE(8,20) 
DO 12 I = Den 
SUMM = 0 


0 
DO 13 J = 1,N 
SUMM = SUMM + AAT(I,J)*SIGMA(J)+GAMA*BBT(I,J) 
13 CONTINUE 
VI(I) = SUMM + COS (TH(T)- ANG) 
I 


Chin Me ue 
WRITE (8, 30) 1 “X(1), ae . ,VI(I),GAMA, SIGMA(I), my 
WRITE (6 EOYs = 8) ie 


12 CONTINUE 
20 FORMAT (//, Be Coe TSP Ol il als 7K, ee YG 'YEL(T) *, oN; 
1 ‘GAMMA! ,354, 'SIGHA 2). 4x, Raton, if? 
30 FORMAT ( 2X Ty, PCTS r 3x ,F8.5,3X, eB. SAV ae 3, SK, Foros 3K,F8. 5) 
SI ~FORMSI 25 or 5) 
RETURN 
END 


PNL(I) 


DOUPWNF OD MOIDMPWNHrH 


eel el pel ae ee ee 


iL; 


ene seselololololelelerelololololelerelelelelelseleloleololelolales 


APPENDIX B 
PROGRAM OUTPUT 


SOURCE AND VORTEX PANEL SOLUTION 
NaCa 23012 
NUMBER. OF PANELS 
ANGLE OF ATTACK 


cE) Aa)’ 
.97500  -0.00350 
.92500  -0.00965 
.85000 -0.01695 
.75000  -0.02580 
.65000 -0.03335 
.55000  -0.03920 
-45000 -0.04325 
.35000 -0.04470 
.27500 -0.04370 
22500 -0.04125 
.17500 -0.03735 
rizsoe -0.03210 
.08750 -0.02765 
.06250 -0.02435 
.03750 -0.01985 
.01875  -0.01470 
.00625 -0.00615 
.00625 0.01335 
.01875 0.03140 
703750 0.04260 
.06250 0.05355 
.08750 0.06115 
.12500 0.06810 
P5008 0.07345 
.22500 0.07550 
.27500 0.07575 
.35000 0.07345 
.45000 0.06775 
.55000 0.05940 
.65000 0.04915 
.75000 0.03720 
.85000 0.02380 
.92500 0.01300 
.97500 0.00460 

EL = 

CD = 

CMLE = 


= 34 
=12.0000 


VEL(1 > 


OR RRR RRRRERERERPMNNE O 


»89112 
»88365 
-88620 
pics) (ew 
-87216 
nooo / 
~82394 
~77104 
L330 
.66069 
-57694 
~46811 
—56096 
e203 92 
OO Ga7 
VATS 
~41005 
gle LlOne)s, 
Besley aie’ 
ere 
17 Biles 
Oe ome 
-90685 
Ri ley S 
sie oy e108) IL 
SD oeos 
~45669 
~ 35494 
se Fleas) 
-20826 
~ 13842 
S010 2 IES) 
HOMME ite 
stele gl 


ieee 
speceeae 
Oui 
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GAMMA 


eee sete eleleroreyerelelerelelelerelelelelololelelolelele’ 


a5 7035 
232035 
moO 5 
22035 
39055 
50 35 
woeUeS 
Cie Oks! she 
soe 
mo 035 
-39025 
sooUSo 
a 2035 
32095 
#39085 
“35055 
a 055 
27055 
wooo 
se Oss 
So 0S 
so 7025 
so 2055 
25.7035 
sey sh0l Sie 
S208 5 
2055 
scUsS 
~o70S5 
Tye Dare 
SSS 
peels 
29055 
OS. 


SIGMA(I) 


PONMNNNNNMNM py pnw pr pon ppp 


408 8 @ @ 6 
MMMM NF Fr 


28747 

31831 
.40999 
.41656 
.42538 
.43692 
.44450 
.47951 
.50195 
54546 
.60630 
.65463 
.64815 
.62380 
.65274 
.58301 
.67224 
.18981 
Echie ts 
.62967 
84536 
.00811 
.23615 
.40177 
.46583 
2.50044 
.58751 
.61244 
.62994 
.65075 
.65258 
.57577 
.59479 


55539 


‘o) 
a 6 
o~ 
tt 
we 


7050 
5 212 ky, 
-21464 
ely ¢ 
e573 
627602 
5s alles 
»40549 
-49107 
56348 
-66714 
1/8087 
yele eal 
Sie ae fe 
S eielhels 
-82974 
226025 
Tolood 
old 
°41720 
~74423 
sodes/ 
-63608 
-0494] 
moO 7 
-42910 
-12194 
5(8) S) eels! 
Aleyaerel 
-45989 
- 29600 
~14958 
-02761 
eeOol!) 


= APPENDIX C 
TABLE I 


TABLE 2 
EFFECT OF GGTRANDXNTIRU ON THE BUGEBE LENGTH 


.6387-7593 






.6567-.6745 .6567-.7093 .6567-.7263 















.6745-.7093 





XTRU=.65| .6567 ‘6 ia -.7420 





-.6920 















.6567-.7263 .6204-.7908 








.6745-7093 


.6745-.7093 6745-.7 | 6387-7429 | .6204-.7752 .6204-.8207 
| NIRU =.68| .6367-.7263 | .6387-.7593 .6387-.7593 .6204-.8060 .6020-.8489 


| NTRU =.69| = .6387-.7429 .6387-.7593 .6204-.7732 .6020-.8350 .3834-.8876 


| NIRU=.70) .6387-.7593 .6204-.7732 .6204-.8060 .6020-.8623 .6204-.8060 


NTRU=.72) = .6204-.7908 .6204-.8060 .6020-.8350 : .3647-.9216 
.9415-.9740 





(turb) 





NTRU =.74 | .6020-.8060 .6020-.8350 .5834-.8623 .5647-.9319 


6020-8350 | .3834-.8623 | .3647-.3752 gl 





NACA 65-213 


Reynolds No. = 240.000 

AOA =(.0 deg 

SE By ls = Begin of transition ( Upper Surface ) 

RIE = Begin of transition ( Lower Surface fixed at 0.4 ) 
Chik = 


Empirical Constant (Gaye) 
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= APPENDIX D 
TABLE II 


TABLE 3 


Pfr eC 1 Ov GG@men ao senC ON THE SHAPE FACTOR(H) 
Porno ZERO SKIN FRICTION 








NTRU =.65 2.745 2.684 2.7024 aeA01 3 





NACA 65-213 


Reynolds No. = 240.000 

AOA = ().() deg 

H =o 0 

NERC = Begin of transition (Upper Surface) 

ae = begin of transition (Lower Surtace fixed at 0.4) 
GG = 


Empirical Constant (Gy) 
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= APPENDIX E 
TABLE Il 


TABLE 4 
EFFECT OF GOGTRAND XATRU ON THE BRAG CObEc I@GIE-yT (CD) 


NTRU =.64 TEsOW2s OLOLS 010029 009876 009813 
WK:.01064 1034 010423 O1LO274 A WO ape, 


NUR = 05 TE:.01022 01013 OTO00zZ2 009900 009785 













NORA aD 01052 OLO416 010301 010193 
XTRU=.66| TE:.01010 010112 O10018 009845 009741 
WK:.01049 01030 010398 010248 010149 

XTRU=-67) Te“ 01010 00926 009785 009816 | 

| WK 1GsS 01030 010316 OLO185 010236 
XTRU=.68| TE:.01G16 009980 009944 009846 009865 
WK:.01055 010383 010338 010246 010312 
| SR Geno TE:.01001 010013 009948 09886 010014 
| WK:.01030 010405 010339 010316 010467 









NTRU =.70 TE: OT 


WK:.01049 


OLOOLS 009900 009964 AU) OME ss 
010403 010334 OLUS82 UNO 29 





TE:.01008 | 00997 009981 010034 010644 
W K:.01047 010332 010380 010481 Ol S2 
ATR =.74 TE:.01004 01005 010017 010489 
W K:01042 010405 _ O10S11 010947 
XTRU =.76 TE:.00999 010024 DIOZ22 
WK 010s? 010434 010672 
NACA 65-213 
Revnolds No. = 240,000 
AOA = 0.0 deg 
"LE = Cd at the trailing edge 
Wiis = Cd at the wake 
OLE Re = Begin of transition (Upper Surface) 
NTRL = Begin of transition (Lower Surface fixed at 0.4) 
GGTR = 


Empirical Constant (Gayep) 
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= AigteeiINiLN F 
TABLE IV 


TABLE 5 
BaRECT OOGIReWND XTRU ON THE ECIFT COEFFICIENT (CL) 








GGTR=10 


NTRU =.64 


| 


| ATRU =.65| es) | .13602 


| 


XTRU =.66 | 


XTRU = | 1596 | 16072 1625 | 1677 : 1800 
| een 1624 | 16464 2074 
=. : / Ls te o> - 
| ERIE . le7l | 1708 1784 | 204 
| 
NTRU =.76 | 1740 | - 1805 1910 : . 
NACA 65-213 
Reynolds No. = 240,000 
AOA =.0 deg 
NOR = Begin of transition (Lpper Surface) 
NRE = Begin of transition (Lower Surtace fixed at 0.4) 
GGTR = 


Empirical Constant (Garey) 
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